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EARTH-MOON ROCKET TRAJECTORIES * 


BY 
LOUIS GOLD! 
ABSTRACT 
This imposing three-body problem is attacked in the approximation of an 
adiabatic limit for which the formalism of securing orbital information is posed. 
An expression is developed for evaluating flight times for rectilinear orbits, but because 
of computational difficulty, an approximation based upon an earlier work is called 
upon. In this manner, lunar flight periods for different launching velocities are 
established. 


1. INTRODUCTION 


With the achievement of rocket propulsion bordering on capability 
for escape flights from the earth, it becomes of immediate practical 
concern to characterize earth—moon trajectories. In the course of such 
a study, it was learned that Egorov (1)? in fact had been investigating 
earth—moon orbits. His approach, however, turned out to be quite 
different from that to be described in this communication ; brute force 
numerical-machine methods were and are being employed by him to 
delineate the various orbits. 

A rigorous treatment of the manifest three body problem is ex- 
tremely difficult. Even the restricted analysis by Egorov (1) requires 
laborious machine calculations. The analysis offered in this report may 
be viewed as a first approximation, being based upon an adiabatic limit 
which considers the earth and moon as essentially stationary during the 
rocket flight. This certainly is an upper limit for launching velocities 
Yo + ©. The chief virtue of this approach, apart from the fact that 
the reduced problem becomes somewhat more tractable, is that it 
affords a means for delineating those regions of flight conditions for 
which more detailed numerical and machine calculations may be con- 


* Presented at the American Rocket Society Semi-Annual Meeting, Los Angeles, Calif., 
June, 1958. 

1 Senior Scientific Executive, Edgerton, Germeshausen & Grier, Boston, Mass. 

? The boldface numbers in parentheses refer to the references appended to this paper. 


(Note—The Franklin Institute is not responsible for the statements and opinions advanced by contributors in 
the JouRNAL.) 
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sidered worthwhile. This is to say that the present study lends itself 
to narrowing down the broad realm of solutions to those of special 
interest, permitting greater concentration of effort where it may count 
the most. 

For further convenience, aerodynamic drag during the early flight 
period is neglected with additional simplification introduced by con- 
fining attention to impulse launching only. Many other factors are 
ignored in this first approach such as the earth’s departure from 
sphericity, etc. (2). 

In the present paper the emphasis will lie in the direction of setting 
forth the formalism and deducing some immediately obtainable findings. 
The more elaborate development will be left to later treatment since the 
problem, even in its reduced form, is still mathematically quite formid- 
able. Thus specific results in this writing will be largely confined to 
properties of rectilinear flights such as flight times and landing velocities. 
In this connection, an upper bound will be drawn in ‘terms of some 
previous work on escape flights from the earth (2). 


y id 


Fic. 1. General spatial configuration for earth-moon flight. M, and M; denote earth 
and moon masses with associated radii Ry and Ry. Rocket position indicated by m; launching 
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2. FORMULATION OF THE ADIABATIC METHOD 


The physical nature of the problem is contained in Fig. 1. The 
perturbing influence of the sun is assumed to be negligible. The sub- 
scripts 1 and 2 refer to quantities in the earth or lunar system. Thus 
a hamiltonian in the adiabatic limit may be written in the form 


GMm 


(1a) 


Lal 


= + 72762") — 
Since the boundary conditions are 
#=0,7; = Rij v= > 0 = (2) 


where vz is the escape velocity, it is sensible to employ the hamiltonian 
(1a) and eliminate 7; as follows. 


dS, ainB, 


Fic. 2. Defining triangle for Eq. 5. 


A; Cos®, 


Referring to Fig. 1, it is seen that 


d= 
0 = 7, cos 6; — £2 COS 82 


whence it is found that 


cos 6; d — r, sin 0, 


fo = 7) tan = 
COS 6; 


cos 6,’ 


Then according to Fig. 2 
r, cos 6; 
Vd? +r? — 2dr; sin 6, 


cos 42 = 


with the outcome 


= (d? 2dr; sin 


| 
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It is to be noted that the treatment is a two-dimensional one in polar 
coordinates; if the earth or moon motion were allowed for, a three- 
dimensional analysis (in spherical coordinates say) would be required. 
Here is an immediate advantage that accrues in the adiabatic limit. 

Dropping the subscripts from r; and 6;, as it will be understood from 
now on that the earth is the reference frame, the Lagrangian has the 
appearance 


GMam 
2 262 


(7) 


Then the equations of motion stemming from 


d 22 29 
dt \ a4; aq: 


turn out to be 


GMm(r—dsin@) _ 
+ ~ ° (8) 
mr*6 + 2mrr6 — cont = (), 


(d? + — 2dr sin 


The solution of (8) is in general difficult. It obviously reduces to the 
standard case in celestial mechanics for M; = 0 and mr*é = const = p. 
Here the angular momentum ? is no longer a constant of the motion. 
For M, = 0, a perturbation solution may suffice; but for the problem 
at hand M,/M, = 80 so that such an attack may give limited accuracy 
in the final results. Elaborate treatment of the general solution is 
deferred for later consideration. 


3. RECTILINEAR EARTH-MOON ROCKET FLIGHTS 


Restricting attention to the rectilinear flight defined by 


t=0, = 3(2n+ = 0 for all (9) 
: the special case for (8) reduces to 
: 1m GMym _ 


Evidently lunar landings correspond to ” even with only the minus sign 


applicable. 
Thus a v(r) relation evolves from the formulation 


[oa = cf (11) 
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which gives rise to 


Now for some value of r a balance between the gravitational pulls 
of earth and moon will prevail. Taking v = 0 at this point means 
setting = vg. Hence 


(13) 


(14) 


which may be re-written in terms of go, the earth’s gravitational constant 
at the surface where r = fy 


= 2g0ro (1 = 2goro. (15) 


It is to be noted that vz is actually slightly less than the ordinary escape 
velocity which of course was to be expected since the moon’s gravi- 
tational field drags the vehicle onward. Defining a quantity ¢ as 
c? = vo? — vz’, Eq. 12 takes the form 


(16) 


where c* is taken always positive for lunar flights. 
Next the minimal launching velocity v, for a rocket to strike the 
moon obtains when r = d — R, with v = 0 from (12): 


To 0 


Then as 
1 1 1 1 


a reasonable approximation for (17) becomes 


R 
(1 ) = gore, 
v + Zoro 


5 
M 1 M 2 ) ( 
M 1 M Be nya a M 1 
with the outcome 
M, M 
Ve 2G ( | 
M 1 M 2 ) : 
as 


6 Louts [J. F. 


where g» is the gravitational constant on the moon’s surface. Com- 
parison with vz in Eq. 15 evidences that 


Vm? + 2g0Fo (19) 
0 


4. CHARACTERIZATION OF DIRECT LUNAR FLIGHTS: 
FLIGHT TIMES AND LANDING VELOCITIES 


Equation 16 lends itself to the calculation of lunar flight times for 
prescribed launching velocities. Thus for rectilinear trajectories 


which yields for ¢(r) 


peel’ Nee (21) 


N: = 2GM.d/c? 


Putting r = d — R, gives the flight time T. 

The evaluation of (21) and the determination of TJ run into the 
complication of nasty elliptic integrals which are not easy to handle 
(see Appendix I). Consequently the rigorous calculation is postponed 
for later consideration and attention is limited to a lower bound solution 
which neglects the moon’s attraction. For this purpose it can be shown 
(Appendix II) that escape flights have the following r(¢) relations: 


V2rogo 

= + (22a) 
2. > 
t = (vo? — Ve (r + — Vro(ro + Co) — colnf(r)} (22d) 


rotcot+ Vro 


(226’) 


Employing the following data 
v* = 7.0 mi./sec., 7% = 3960 mi., r & 240,000 mi., 


a rocket dispatched with v) = v* ought to take about 


where 
Ve? r+o+ vr 
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T= (24)#/2 & 10® sec. = 1.8 X sec. 59 hr. 


=3 7 (3960)12 


Cole (3) has made such an estimate and gets = 49 hr. The present 
method can next be utilized for flights where v» > v*; a reduced form 
of (226) is adequate here 


(22c) 


This gives rise to further data in Table I. 


TABLE I. 


T, hr. 50 5.5 1 
1 2 10 


The 50-hr. flight associates as noted to a zero landing velocity. The 
landing velocities for the more rapid transits may be calculated from (2) 


v*? — (vo? — vt? — yt 


(Wo? — + 00) (Vo? — — v) 
— — 09) — v*? + 


(23) 


Z= 


APPENDIX I 
To make use of Byrd & Friedman Tables (4) the integral in (21) must be put into the form 


— — d) 


Then for all constants real there are four cases 


a>b>c2r>d, (4, p. 103) 
a>b>c>re2d, (4, p. 107) 
a2zr>b>c>d, (4, p. 120) 
a>r2zb>c>d, (4, p. 124). 


Additional cases arise when any of constants complex. 
A typical illustration has 


(¢ — c)(t — d) 1 dn?u 


a—b (a — c)(t — d) 


g= 2(a d)-"2, 


ve 
where 2 


8 Louis GoLp 
Further evaluation encounters for this special case 


dn*u 1 


atsnucnudnu 

+ (2a? — at — k*) + — | (I.3a) 

with 

(a — c)(b — d) 


(1.38) 


Thus it turns out that more effective ways of computing ¢ are to be found. Several schemes 
for doing so will be dealt with in a later work. 
APPENDIX II 


A calculation paralleling the one made here may be found elsewhere (5); the method 
given below is more explicit and does not contain peculiar, restrictive boundary conditions. 
Starting with 


it is easy to find 


(II.2) 


Three cases must be considered : 
1. v9 = v* (critical lunar flight) 
1 


1 2 


which yields (22a). 
2. v9 < v* (earth bound flights) 


r 
whence 


t = (v*? — rola —ro) — Vriq—r) +4 — sin 
ve? 
q=To (i 
3. v > v* (lunar fights) 
The result (226) follows directly from (1.4) in this situation. 
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EFFECTS OF UNIPOLAR AIR IONS ON MICROORGANISMS 
AND ON EVAPORATION * 


BY 
A. P. KRUEGER,' W. WESLEY HICKS,’ AND J. C. BECKETT * 


SUMMARY 


Bacteria in droplets of distilled water were exposed to unipolar atmospheric ions 
produced by radioisotopes. Experiments using micro-chemical equipment and bac- 
teriological techniques were carried out both in clean and contaminated air. Lethal 
results in clean air were attributed primarily to the action of ions on staphylococci. 
Survival in contaminated air was attributed primarily to the action of contaminants 
in shielding the staphylococci from (high mobility) ions. Evaporation rates were 
found to differ according to the presence, or absence, of air ions and atmospheric 
pollutants. 


Although it is well recognized that certain properties of the ambient 
atmosphere, for example, humidity, greatly influence the fate of air- 
borne microorganisms, the effect of atmospheric electricity has been 
little more than suspected. At times, conditions would seem to favor 
the existence of such action for air ions can be present in the atmosphere 
in considerable concentrations and in various sizes, depending upon pre- 
vailing weather conditions, air pollution and ionizing sources. While 
dried bacteria carried in the air are relatively resistant to most climatic 
factors (heat, for example), moist cells are considerably more susceptible 
to both physical and chemical forces. In the present work, micrococci 
suspended in small drops of distilled water were exposed under con- 
trolled conditions to unipolar air ions generated by alpha or beta radi- 
ation derived from radioisotopes. 


THE BACTERIOLOGICAL EXPERIMENT 


During the past two years, Krueger, Smith and Go (1)* working in 
Berkeley, California have conducted bacteriological experiments to de- 
tect air ion effects on bacteria. Preliminary work failed to reveal any 
influence of air ions on such activities as the swarming of Proteus vul- 
garis, the production of hemolysin by streptococci, or the growth rates 
of various organisms in nutrient media. However, it was found that 
significant changes in the survival time of staphylococci were produced 
if the experimental design included the following points: 

* Presented before the First International Bioclimatological Congress of the International 
Society of Bioclimatology and Biometeorology, held in Vienna in September, 1957. 

1 Department of Bacteriology, University of California, Berkeley, Calif. 

2 President, Wesix, San Francisco, Calif. 

3 Chief Engineer, Wesix, San Francisco, Calif. 

‘ The boldface numbers in parentheses refer to the references appended to this paper. 
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1. Exposure of cells in droplets sufficiently small to provide a high 
ratio of surface area to volume. 

2. Suspension of washed organisms in distilled water rather than 
nutrient fluids. 

3. Adaptation of microchemical equipment to bacteriological tech- 
niques for handling bacteria in small drops and for quantitative 
estimation of viable cells in fluid aliquots varying from 4 to 100 
lambdas (lambda = 10~ liters). 
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Staphylococci were suspended in distilled water and droplets con- 
taining 4 to 100 lambdas of the suspension were placed in a glazed ceramic 
crucible. Ions were produced in the air above the crucible by radio- 
active means, and the air ions were directed toward the drops in the 
crucible by an electrostatic field. After exposure to ionized air for 
various intervals of time, the remaining living organisms were counted 
by means of replicate plate counts. Results obtained under conditions 
of clean air were as follows: 


1. Positive and negative ions significantly increased the rate of cell 
death (Fig. 1). 

2. No important changes of pH or Eh occurred. 

3. The rate of evaporation from the droplets was increased slightly 
(Fig. 2). 


The polonium and tritium sources were used, but by chance the 
210Po experiments were done primarily during the period of greatest air 
pollution (August to November) and the *H experiments at a time when 
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the air was relatively clean. Differences observed were attributed to 
the quality of the ambient atmosphere and not the ion source. 

In the tritium experiments performed in the absence of air pollutants, 
the death rate of exposed staphylococci was measurably increased and 
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so was the rate of evaporation. The lethal effect was found to be pro- 
portional to the numbers of air ions reaching the surface of the droplet. 
At 16-cm. distance only 5 X 10‘ air ions/mm.?/sec. arrive at the surface 
and this is too small a quantity to induce any measurable biological 
changes. 

The observed increase in rate of evaporation appears to be a reason- 
able consequence of air ion action, for the motion of air ions in a charged 
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field actually constitutes a minute electrostatic wind and one would 
anticipate an increase in the kinetic energy of water molecules in the 
suface layer of the droplets. The tritium experiments further suggest 
that in a clean atmosphere air ions can act directly on cells. Thus, if 
the total volume of the droplet was kept constant through the addition 
of distilled water at intervals to prevent effects due to evaporation and 
the droplet was stirred to facilitate contact of bacteria with air ions, 
pronounced lethal effects were obtained while unstirred controls exhib- 
ited no changes from the normal rate of cell death (Fig. 3). 
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FLUID ADDED AT TIMES INDICATED BY ARROWS 
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The general features of the “Po series apply to 100 droplets of a 
distilled water suspension of staphylococci in smog-polluted air and may 
be summarized as follows: 

1. Positive and negative ions delayed the fall in pH and accompany- 

ing rise in Eh customarily observed in smog-laden air (Figs. 4 
and 5). 
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2. Both ions diminished the rate of water loss from exposed droplets 
(Fig. 6). 

3. The result of 1 and 2 was a protective action on suspended bac- 
teria; that is, the death rate of staphylococci was diminished 
(Fig. 7). 

4. In air saturated with water vapor the air ions exerted no meas- 
urable effect. 

It seems probable that the protective effect observed in the polonium 


experiments was due to the absorption of positive or negative charges 
by smog particles, with the result that they became mutually repellant. 
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Accordingly, smaller numbers of smog particles per unit of time arrived 
at the water droplet—air interface and the fall in pH usually occurring 
as the smog constituents dissolved in the droplet was delayed. Also 
operating to protect the suspended organisms was the lower rate of 
evaporation that prevailed in preparations exposed to air ions. It is 
known that large (Langevin) ions form continuously as small air ions 
act on certain smog constituents. Such large ions can serve as centers 
of condensation for water vapor and the water thus contributed to the 
droplet could replace to some degree the fluid lost by evaporation. If 
this mechanism were to function without inducing a change in pH it 
must be assumed either that the smog components involved in the for- 
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mation of Langevin ions are not acidic in character or that the numbers 
serving as centers of condensation are too few to affect the pH of the 


droplet. 
GENERATION OF IONS FOR BACTERIOLOGICAL STUDIES 


Of the several sources available for ion generation in the laboratory, 
radioisotopes have been judged the most suitable, since they function 
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without producing by-products such asozone. Worden (2), for example, 
has successfully used a sealed source of polonium 210 in experiments on 
tissue cultures. Radiation from polonium is limited to alpha particles 
travelling a maximum distance of 4 cm. and extremely weak gamma 
radiation not detectable with ordinary instruments. 

Unipolar ions are produced by separating positive and negative ions 
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in an electrostatic field (3). The single-polarity ions are directed to the 
target area using the same electrostatic field, Fig. 8. By varying the 
distance from the source material to the target and the potential be- 
tween source and target, the number of unipolar ions reaching the target 
can be accurately controlled. 

Subsequent to the development of the polonium-type ionizing tube, 
tritium foil (titanium containing tritium gas) became available and was 
also used in these experiments. Tritium has the advantage over polo- 
nium of a much longer half-life—thirteen years versus 140 days. 
Furthermore, it is possible to concentrate tritium in a sealed foil source 
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Fic. 8. Arrangement for exposure of bacteria to unipolar ions. 


to a much greater degree and still retain a high factor of safety from 
leakage of the radioactive material. For example, polonium ionizing 
tubes normally use 6 sq. cm. of foil and contain one-half millicurie of 
polonium. Tritium units are constructed using 0.4 sq. cm. of foil and 
contain 50 millicuries of tritium. 

Tritium is a soft beta emitter, 0.015 mev, with no alpha emission 
and an insignificant amount of secondary radiation. Typical values of 
Bremsstrahlung (gamma) readings from a tritium ion generator are 
5MR/hr. at 6.3 cm., 2.3 MR/hr. at 11.4 cm., and 0.8MR/hr. at 16.5 
cm. The small size of the tritium foil itself permits the use of a smaller 
structure for ion separation and this is of real advantage in studies where 
space is limited (4). Low energy beta radiation is not as efficient for 
ionizing air as alpha particles but the greater ease with which tritium 
can be concentrated in sealed form permits manufacture of more effec- 
tive tritium ion generators where high density is desired. The relation- 
ship of ion density to distance and field potential is shown in Fig. 9. 
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ION MEASUREMENT 


Measurement of atmospheric ions can be made in various ways (5), 
but for this work a simple target collector was selected as being the 
most informative. The bacteria are suspended in small droplets in a 
glazed ceramic crucible and ions are directed toward the crucible as 
charged particles toward a target. 

A metal disk, 25 mm. in diameter, forms the target of the measuring 
apparatus, Fig. 10. This disk is surrounded by a grounded wire shield 
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Fic. 9. Unipolar ion density (tritium source 50 MC) as a function 
of distance and potential. 


so that the target will not distort the electrostatic field. The disk is 
attached to a well insulated and shielded conductor for connection to a 
micro-microammeter. Ions reaching the target area are measured in 
terms of ions/mm.?/sec. as determined from the equation N = I/gA 
where: 


N = number of ions/mm.?/sec. 

I = ion current in amperes 

q = charge of one ion—1.6 X 10—" coulomb 
A = area of probe in mm.’ 
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Results using the target collector in an unconfined space are illus- 
trated in the chart, Fig. 9. It is chiefly high mobility ions that reach 
a target collector, although this device measures the charge of all ions 
reaching the target, regardless of size. 

These experiments were designed to detect possible lethal effects of 
air ions, and these densities are substantially higher than found in 
ordinary outdoor air. With about 200 microcuries of polonium foil the 
air ion density inside the exposure chamber was observed to be 1 « 10° 
ions/mm.?/sec. at a working distance of 7.5 cm. and 300 volts. In 
similar experiments using a 50-millicurie tritium source, the correspond- 
ing figures of 4 X 10° ions/mm.*/sec. at 4 cm. were measured. 
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Fic. 10. Ion current probe. 


Appropriate steps were taken to avoid two possible sources of error 
in this type of experimentation: namely, direct action of radiation on 
the cells, and distortion of viability counts by bacterial agglutination. 
The distances between the polonium or tritium foils and the test drop- 
lets are much greater than the calculated maximal ranges of alpha or 
beta radiation. Bremsstrahlung or secondary radiation is extremely 
low from tritium. Likewise, gamma radiation from polonium is barely 
detectable. As a precaution, control experiments were run with each 
radioisotope present but without the electrostatic field to effect ion 
separation. Similarly, experiments were made with the electrostatic 
field and no radioisotope to insure that neither by itself had any signifi- 
cant effect on the organism. Cell clumping is not observed by micro- 
scopic examination and, further, the time periods required to induce 
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lethal effects in light cell suspensions are less than the minimal times 
for agglutination predicted by the von Smoluchowski equation (6) : 


t = 
where 
i = time in seconds 
r = radius of particle 
Vo = number of particles/ml. 


and 
D = 


is the Einstein equation (7) for the diffusion coefficient. 

An independent study of air ion effects has been conducted recently 
by Dr. Robert Fuerst and R. J. Bell, Jr., at the University of Texas, 
M. D. Anderson Hospital and Tumor Institute. Using several strains 
of Neurospora crassa they found that the germination of spores was 
inhibited by exposure to either positive or negative airions. As in the 
experiments on staphylococci reported in this paper, fairly prolonged 
periods of exposure were required to induce the effect. 


CONCLUSION 


First, in clean air high mobility ions derived from radioisotopes have 
a lethal effect on certain bacteria suspended in droplets. Second, in 
clean air high mobility ions increase the rate of evaporation of droplets. 
If the air contains pollutants, the reverse occurs in each case. 
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Floating Zone Refining of Thermal 
Sensitive Intermetallic Compounds.— 
An improved technique for growing 
single crystals of binary semiconduc- 
tors which decompose on melting has 
been described by J. M. Whelan of 
Bell Telephone Laboratories. The 
basic experimental work was _ per- 
formed on gallium arsenide, but the 
method should be applicable to a va- 
riety of compounds which are ther- 
mally unstable at their melting points. 

According to the paper presented 
by Dr. Whelan to the 133rd National 
Meeting of the American Chemical 
Society in San Francisco, the floating 
zone method appears to be superior 
to other methods for growing GaAs 
single crystals. Composition of the 
liquid phase at the melting point is 
strongly dependent on the partial 
pressure of arsenic. This is most 
easily controlled by using a sealed 
system containing excess arsenic and 
regulating its minimum temperature. 
The rate of growth is easily controlled 
and spurious nucleation is greatly re- 
duced by the thermal symmetry of the 
freezing interface. 

Other advantages of the method are 
the reproducibility, regularity of im- 
purity distributions, and the minimi- 
zation of apparatus contamination. 
The relatively large surface-to-volume 
ratio favors the approach to equilib- 
rium between the liquid and vapor 
phases. In addition, the small liquid 
volume minimizes temperature varia- 
tions in the melt, and the consequent 
concentration gradients of the princi- 
pal components. 

The technique should be most use- 
ful with binary compounds in which 
only one of the component elements 
has a considerable vapor pressure at 
its melting point. The compound 
must have a high enough electrical 
conductivity to allow heating by radio 
frequency induction. The _ surface 
tension and density of the molten ma- 


terial must also be such as to support 
a molten zone during the process. 

In the basic floating zone refining 
technique, a rod is supported verti- 
cally. A heat source, for example an 
induction coil operated at radio fre- 
quencies, is moved relative to the rod, 
melting a liquid zone as it moves. 
Surface tension supports the liquid 
zone. Usually by this method, a 
single crystal can be grown and puri- 
fication achieved in the zone refining. 


Supercharged Hydraulic Oil Sys- 
tem Eliminates Surging.—The Yale 
& Towne Manufacturing Company 
has announced the perfection of a 
supercharged hydraulic oil system for 
lift trucks to eliminate surging in the 
hydraulic tank and the danger of get- 
ting air bubbles into the lines and 
pump. This new development is be- 
ing made a standard integrated de- 
sign feature in all Yale rider-type gas, 
LP-Gas and electric powered lift 
trucks. 

The presence of air in excessive 
amounts in a lift truck hydraulic sys- 
tem can lower hoist speeds, cause 
hoist cavitation and, in extreme cases, 
lead to pump failure. In the new 
Yale system, the free flow hydraulic 
return line brings the hydraulic oil 
into the tank at a midpoint on the in- 
board side. A metal pipe conveys the 
oil straight down into one end of a 
channel-shaped passage welded to the 
bottom of the tank. The suction line 
leading up out of the tank into the 
hydraulic system is located at the 
other end of this passage. Action 
of the hydraulic pump pulls the oil 
through the passage and into the suc- 
tion line in a smooth, steady flow. 
A series of apertures along the sides 
of the channel passage bleed off ex- 
cessive oil and pressure smoothly, pre- 
venting a bubbling action in the oil 
remaining in the tank. 
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PLASTICS IN THE SPACE AGE * 


BY 
JOHN H. LUX? 


On October 4, 1957, our radios told us that the Russians had ushered 
in the Space Age with the first satellite. We all knew it was coming, 
but we didn’t expect it so soon, nor did we expect to find ourselves 
behind in the ‘Space Race.” 

This has required all of us to evaluate our own future contributions 
so that we can do everything possible to obtain a lead in the near future. 
We feel that the plastic industry will make major contributions to the 
“Space Age.” 

Because of the widespread use of plastics in every phase of our life, 
many people have suggested we were in the Plastic Age. We have 
plastics in clothing, automobiles, homes, appliances and food packaging. 
The word “plastics” brings many different images to peoples’ minds, 
yet to 99 people in a 100, plastics means plates, drinking glasses, tooth- 
brush handles, automobile trim, kitchen work surfaces, toys, or swim- 
ming pools. These products have emerged from the laboratories in 
recent years and are providing us with new and useful items, many of 
which we never enjoyed before. Plastics in the coming “Space Age”’ 
will multiply this list and add many more items and uses, some of which 
are beyond the scope of our present thinking. 

Let’s look at the demands of the future. First, there will be a 
demand for versatile materials of construction for making huge struc- 
tures which are lighter in weight, unaffected by the elements, easy to 
handle, transparent or translucent. Second, there will be a demand 
for materials to operate at temperatures lower than those now generally 
utilized. Third, there will be a demand for materials to withstand the 
high temperatures—25,000° F. and perhaps higher. Fourth, there will 
be a demand for materials which will be expected to stand the high 
erosion of re-entry and the stress and strains of speeds of high Mach 
numbers. Fifth, products will be needed for electronic applications 
with a wide variety of electrical, magnetic and transistor properties. 
Sixth, there will be a demand for materials which will be able to with- 
stand the nuclear radiation, cosmic rays and other problems encountered 
in space travel. Seventh, we will need plastics which will provide large 
amounts of energy through combustion or other chemical reactions for 
propulsion. 

* Adapted from a lecture delivered at the Stated Meeting of The Franklin Institute, 


Wednesday, April 18, 1958. 
1 President, Haveg Industries, Inc., Wilmington, Delaware. 
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Many of these properties are now available to a wide degree in plastic 
materials but we look forward to vast improvements in the future. 
Before looking at plastics, let’s look at why the other materials are 
unsatisfactory under the above conditions. In the high temperature 
and high velocity field (areas above 2000° F. and speeds over Mach 1) 
most people felt that the answer to high temperature materials of con- 
struction had the best possibility of coming from such structures as 
metals, metal carbides, graphite or ceramics. With respect to metals, 
the thermal conductivity is so high that the metals very quickly reach 
the temperature at which they lose strength and even melt or flow. 
Sometimes this happens as quickly as a few parts of a second, particu- 
larly above 3000° F., so that above this temperature metals are not 


High temperature heat exchanger system. 


Fic. 1. Nose cone. Fic. 2. 


practical. In the case of carbides or graphite, these materials at similar 
high temperatures are subject to oxidation in an oxidizing atmosphere and 
are also subject to reduction in a reducing atmosphere. Carbides and 
graphite also have a very high thermal conductivity and both of these 
materials are subject to serious limitations as far as fabrication tech- 
niques are concerned. No large structures can be made without great 
difficulty and there are limitations as to the strength of even smaller 
structures. In the case of ceramics, these materials are subject to shock 
so that in going from room temperatures to 5000° F. in a matter of a 
second or even less, these materials are subject to cracking. 

In contrast to this, plastic materials and their derivatives have 
resistance to oxidation and reduction, good thermal insulating qualities, 
ability to withstand shock, ease of production by a wide variety of 
techniques, so that they can be fabricated from } ounce up to 50,000 
pounds and will operate permanently at temperatures up to 550° F., 
they will operate for thousands of hours in the 600-700° F. range and 
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at lesser times up to 18,000 or 20,000° F. As a result, in the high tem- 
perature field, it has been possible to make such materials of construction 
into blast tubes, re-entry nose cones for missiles, columns, 6000° F. heat 
insulation, heat exchangers, nozzles, aerodynamic fins, missile jet fins, 
jet stabilizing rods and many other thermal structural components for 
aircraft, rockets, missiles and atomic energy processes. (See Figs. 
1, 2 and 3.) 

In the nuclear field, plastic thermosetting materials have been tested 
at 10° roentgens for periods of several months and an examination of 
test samples showed their strength doubled over a period of one year. 
These are now used in equipment for uranium refining. 

In the field of polyethylene resins, it has been found that radiation 
causes a change in structure which will allow the use of polyethylene to 
be upgraded from around 130° F. to temperatures above the boiling 
point of water so they can be sterilized, thus opening up broad new uses. 


Fic. 3. Large nozzle. 


In the low temperature field, the use of Teflon and Silicone rubber 
products enables us to produce materials which will operate at —40° 
and —60° F. and still retain the ability to withstand shock and not be 
subjected to fatigue as would metals. The silicone polymer oils provide 
lubrication as low as —60° F. 

The reinforced polyester plastics were chosen for the large structures 
to house the DEW line, early warning radar system for use near the 
Arctic Circle, because of their ease of assembly, light weight, superior 
insulating and electrical properties. Many future structures, homes, 
factories, automobiles, aircraft, and auditoriums will soon burst forth 
in this versatile combination of plastics and fiberglass. 

In the electronic field, Teflon coated wire makes it possible to en- 
gineer missile systems which will operate at high temperatures. It also 
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enables computers and other delicate equipment to be built to operate 
at temperatures which would degrade most insulations. 

These things have not been necessary in radios, television and other 
products because it is easy to call a repair man for a few dollars and 
replace any defective system. However, as one man has put it: ‘‘At 
20,000 ft. or 20 miles up there are no repair men.” Teflon-coated wire 
can be welded because in decomposition it does not form carbon particles 
which would set up conducting points. In the electronic area the use 
of plastic radomes, electrical terminal insulators, various complex elec- 
tronic equipment, coaxial cabling and many other items are fabricated 
from plastics. 

The guidance and control system of the missile represents about 50 
per cent of its cost (but only 10 per cent of the weight), and it is here 
that plastics really offer advantages. Plastics are used to make the 
printed circuits, in eliminating metal parts, for potting and impregna- 
tion, and protection of wire and electronic components (see Fig. 4). In 


Fic. 4. Antenna covers. 


the mechanical missile components area there is a wide variety of cams, 
gears, levers and structural assemblies. In the manufacturing pro- 
cesses, plastics are ideal for tooling, since design changes can be made 
more quickly and at lower cost than with steel tooling. 

Plastic foams, fillers and honeycomb are used in sandwich type con- 
struction to provide stronger units. Many metal units, such as the 
Martin Matador (B-61) are assembled using plastic adhesives. The 
number of uses of plastic coatings, sealers and liners on other materials 
will be broadened considerably as new materials are developed. 

Also, another feature of the “Space Age”’ will be more difficult cor- 
rosion problems in the manufacture and handling of materials in the 
future. Plastics, unlike metals, do not corrode, so that they can be 
made to handle such violent materials as hydrofluoric, hydrochloric and 
sulfuric acids; strong caustics, chlorine dioxide, aqua regia and chlorine 
gas. The preparation of uranium, columbium, molybdenum, tungsten 
and tantalum for “Space Age’’ use, uses large amounts of plastic 
equipment. 
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The development of plastics for the uses of the past and improve- 
ments for the future comes about only through the expenditure of a 
considerable amount of money and technical effort. This does not 
mean, however, that only large firms can play a significant role in the 
plastics industry. 

My own firm, Haveg Industries, Inc., is a small company, working 
with plastics since the early 1930’s. We have concentrated on solving 
over one million corrosion problems in the last 25 years, and we have a 
number of firsts in the industry to our credit: (1) the first plastic pipe 
(some 25 years ago); (2) the first plastic fittings and valves; (3) the 
first chemically welded field systems; (4) the first plastic tanks (up to 
35,000 gallons); and (5) the world’s largest plastic structure, a fume 
stack weighing 42,000 pounds. This stack, shown in Fig. 5, was field 


Fic. 5. 200-ft. high temperature stack. 


constructed to form a unit 200 feet tall, 5 feet in diameter, designed to 
withstand a 125-mile wind load. It has been service for several years 
replacing stainless steel stacks lasting only six months, and at a money 
saving. 

Using over thirty different plastic combinations, we manufacture 
thousands of items such as tanks, hoods, fume ducts, heat exchangers, 
cooling towers, reaction vessels, filters, steam jet ejectors, drying trays, 
electroplating tanks, bifurcators, pumps, rolls, corrosion resistant 
cements and tank trucks. 

It was this background in plastics that caused Haveg to initiate work 
in the high temperature field many years ago. We felt that plastics 
were usable at temperatures of 5000° F., for short periods of time—and 
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we proved it in actual rocket-firing tests. Our proprietary materials 
known as Rocketon, Missileon, Mercuron, Planeton, Galaxon and 
Asteron for high temperature applications find use in blast tubes, nose 
cones, nozzles, jet vanes, tail cones, insulators and adaptors (see Fig. 6). 
In addition to these fabricated materials, we also produce Siltemp glass 
which is available in both a sheet form and in glass rovings. This 
material is stable at 1800—2000° F., for use as an ingredient in many 
military projects and for solving many high temperature problems. 
Materials can be fabricated for short term use and insulation up to 
6000° F. The latter products are usable in areas needing high resistance 
to erosion. 

One of the most outstanding contributions of plastics to our future 
in the ‘Space Age”’ is that of furnishing the energy for our solid propel- 
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lant type missiles. This field, pioneered by Thiokol Chemical Corpo- 
ration, has changed rocketry from a “‘hit and miss” proposition to a 
field of high reliability. The use of plastics and an oxidizer such as 
ammonium perchlorate provides good combustion properties, ease of 
handling, high strength, and ability to work between —40° and +160° F. 
Among the plastics which are used are polysulfides, polyethylene, epoxy, 
vinyl, polyurethanes and various synthetic rubbers. 

It is obvious that plastics are able to provide the entire missile sys- 
tem with the exception of the wiring and oxidizer; and plastics have 
been used in the manufacture of both of these. 

We have shown you some of the areas of plastics in our future, a 
future in which even a small organization can play a very significant 
part and in doing so, can multiply its business from that of a small 
business to a recognized growth industry in the short period of two to 
three years. 
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THE LAGRANGE EQUATIONS IN ELECTRICAL NETWORKS 


BY 
FREDERICK L. RYDER ' 


SYNOPSIS 


The dynamic equations of Lagrange, which are known to apply for certain 
lumped-element electrical networks, are derived for such networks by a method which 
results in expressions in which the usual differentiations need not be carried out. 
These expressions, one suitable for mesh analysis and the other for nodal analysis, 
are shown to reduce to equations previously obtained by energy considerations in the 
steady-state case, and, like those equations, are particularly useful when transformers 
are present. An example is given to illustrate the method of application, as well as to 
show the relative advantages of the mesh and nodal forms of the expressions. 


INTRODUCTION 


The familiar Lagrange equations for the analysis of dynamic systems 
are widely used in mechanical problems, where they yield the following 
advantages : 


1. The derivation of the governing equations of the system is rela- 
tively simple because it is formalized once the energy expressions are 
known. A further, though secondary, aid to simplicity lies in the rela- 
tive ease of keeping track of polarities, since kinetic and elastic energies 
are always positive. 

2. The solution of the governing equations is relatively compact in 
that the effects of elements of the mechanical system which absorb no 
net energy, such as ideal levers, are effectively eliminated from the 
analysis at an early stage, so that during most of the work the number 
of parameters is less than in other methods. 


It has been shown (1)? that the parameters in certain lumped-ele- 
ment electrical networks can be arranged so that Lagrange’s equations 
are satisfied. However, since such networks lack the multidimension- 
ality of mechanical systems, the Lagrange equations offer no advantage 
over conventional analysis by Kirchhoff’s laws—with one exception, 
namely, when transformers are present. In that case each transformer 
may be represented by a combination of an ideal transformer and suit- 
able circuit elements by a number of different conventional methods, 
after which the Lagrange analysis is found easier than the Kirchhoff 
analysis because the effects of the ideal transformers, which absorb no 
net energy, may be eliminated early in the analysis. 

The electrical form of Lagrange’s equations has a further use in 


1 Member, Scientific Research Staff, Republic Aviation Corporation, Farmingdale, N. Y. 
2 The boldface numbers in parentheses refer to the references appended to this paper. 
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connection with certain electrical-structural analogies (2) which are 
based on the similarity between corresponding electrical and structural 
expressions for the energy terms in these equations. 

In the present paper we shall arrive at electrical forms for the 
Lagrange equations by a method quite different than that used by 
Jeans, namely by the use of Kirchhoff's laws alone, first by a mesh type 
of analysis and then by a nodal type. We shall also develop expressions 
which, while they are equivalent to those of Lagrange, are easier to 
apply in that the numerous differentiations required in the Lagrange 
method have already been performed. 


MESH ANALYSIS 


Let Fig. 1 represent a typical branch jk in an electrical network, 
with constant values of series inductance, resistance and capacitance as 


j 


Fic. 1. Fic. 2. 


indicated. A series voltage difference which rises by the prescribed 
instantaneous value v, in the direction 7 to k is shown, as well as a 
branch charge of instantaneous value ¢;, defined as the time-integral 
of current flowing in the reference direction. For definiteness we shall 
stipulate that gj, and vj, denote the respective values of charge and 
voltage over. and above the initial values when time ¢ equals zero. 
Since charges and voltages are related by linear algebraic or differential 
equations, the principle of superposition assures that initial and time- 
dependent values for any charge or voltage may be added to find the 
total instantaneous value. 

Figure 2 shows a simple network containing several branches, with 
Z denoting the combination of L, R and C for brevity. The generator 
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in branch 1-3 is a source of prescribed voltage of the indicated instan- 
taneous value, the charge g:; being unknown. An additional type of 
element has been added in branch 2—4 in the form of a generator of 
prescribed charge of instantaneous value g2,, the voltage rise being 
dependent on the remainder of the network. (Such a generator can be 
simply made up by placing a suitable voltage source in series with a 
resistance whose voltage drop is large compared with the voltage differ- 
ences across the other circuit elements.) Note that nodes 2 and 3 have 
been given separate designations even though they are at the same 
voltage, so that the various circuit elements may be referred to un- 
ambiguously by means of suitable subscripts; this can always be done 
no matter how complicated the network. Also note that the network 
is complete within itself; if we wished to enlarge it, say by inserting a 
current at node 1 and withdrawing it at 2, the complete additional 
branch would be shown as part of the enlarged network. 

By the conventional mesh technique, with mesh charges substituted 
for the usual mesh currents, all branch charges can be evaluated in 
terms of the prescribed charge plus an additional arbitrarily selected 
mesh charge. In Fig. 2 the additional mesh charge is shown as 41s, 
which together with g2, gives two mesh charges as shown, with the 
result that 

= 
diz = Yea — Q13- 


The unknown mesh charge gi; may be considered redundant in that 
it cannot be determined by Kirchhoff’s first law (of current-continuity) 
alone, in the same sense that structural forces and moments are con- 
sidered redundant when they cannot be determined by equilibrium 
considerations alone. 

These ideas can be extended to more complicated networks, where 
in general any branch charge ¢;, can be expressed as follows: 


= + b jx2Q2 D ms (1) 


where m is the number of meshes, and the g’s on the right contain all 
the prescribed charges and as many redundant charges as are necessary. 
The d’s must be either + 1, — 1 or zero, and can be determined either 
by the conventional mesh technique or by applying Kirchhoff’s first law 
to each node in turn. 

The total voltage drop in the reference direction in the typical 
branch of Fig. 1 is as follows: 


+ Radix + — 0% = — (2) 


where the dot represents differentiation with respect to time, and Z» 
is an appropriate operator. Alternatively, v;, could denote the voltage 
rise across a charge generator. Summing to zero the corresponding 
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voltage drops, expressed in terms of only prescribed and redundant 


charges, around the meshes indicated in Fig. 2, we get 

Z13(G13) Z12(Ges = Gis) = 0 

+ Z 41 (Ges) + is) = 0, 
where v2, is the unknown instantaneous voltage rise across the charge 
generator in the reference direction. These equations could have been 

expressed as follows: 

ix (Gin) — Vx] = 0 s = (3) 
where the summation is taken over all branches for each value of s, 
which represents each of the m prescribed and redundant currents in 
turn. For any particular value of s the }’s assure that only those 
branches which form part of the mesh are counted, and also that the 
direction in which each branch voltage drop is measured is suitably 
taken into account. 


Equation 3 expresses the mesh summations of voltage drop in any 
network, no matter how complicated, made up of the elements shown 
in Figs. 1 and 2. The resulting m equations are just sufficient to permit 
evaluation of the redundant charges as well as the unknown voltage 
rises associated with the charge generators. Since the 0’s will equal 
zero when s denotes a redundant charge and 7k a branch containing a 
prescribed charge, the redundants can always be evaluated inde- 
pendently of the voltage rises of the charge generators. 

If transformers are present, each of them should first be represented 
by a combination of an ideal transformer with suitable circuit elements 
in any one of the numerous conventional ways, after which the current- 
relationship (which is the same as the charge-relationship if the initial 
values of charge are ignored) imposed by each ideal transformer on the 
currents through its various windings should be taken into account in 
writing Eq.1. In the particular case of the ideal transformer illustrated 
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in Fig. 3, where a common flux links three windings having the respective 
numbers of turns N,, N2, NV; and carrying the respective currents 4,, 72, 
43, the fact that the net magnetomotive force acting on the core must 
be zero leads to the following relationship : 


+ + Nai; = 0, (4) 


where the polarity of each of the N’s is such that the associated current, 
when positive, produces flux in a common prescribed direction. The 
effect of the turns ratios on the current-relationships is shown in Fig. 3, 
where 7; is expressed as a function of 7; and 7:2. 

Equation 1 cannot now be developed by straightforward application 
of the mesh symbolism, since in the latter the various transformer wind- 
ing currents would generally be assumed independent because they 
would fall in different meshes, whereas they are actually related by 
expressions such as Eq. 4. However, Eq. 1 can always be developed 
by considerations of current-continuity at each node. In general this 
will affect the 6’s in Eq. 1 so that they may have any constant values, 
depending on the turns ratios. 

To show that Eqs. 3 are valid when ideal transformers are present 
we shall write it for the case s = 1, assuming that the transformer of 
Fig. 3 is the only one appearing in the network. We get 


Db jail Z se (Qin) — Vie] + 1(— — (— es) = 0, 


where the summation term is taken over the entire network except for 
the transformer windings, and the e’s represent transformer voltage 
rises (measured over and above their values at ¢ = 0) as defined in the 
figure. Since induced voltage is proportional to the number of turns, 
we may write 

= = e3/Ns (5) 


so that the last two terms on the left in the preceding equation sum to 
zero. This argument can be extended to any configuration of ideal 
transformers, so that in general the voltage rises of ideal transformers 
can be omitted in applying Eq. 3. This may yield a considerable saving 
over conventional methods of analysis, where the e’s must be carried 
along as unknowns even if only the charges (or currents) are to be de- 
termined. In the present case the e’s are readily calculated, if desired, 
by voltage-compatibility considerations after the charges have been 
determined. 

In the special case where the prescribed charges and voltage rises 
are steady-state sinusoids at a given frequency, it is conventional to 
deal with current vectors I rather than the corresponding instantaneous 
values of charge, and to denote vectorial voltage rises by V and complex 
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impedances by Z. Then Eqs. 1 and 3 may be replaced by the following : 
I. + D + + D (6) 
TZ Vin} = 0 = 1, 2. hie (7) 


The first of these is identical to Eq. 18 in a previous paper (3) dealing 
with energy considerations in the steady state; the second takes the 
place of both Eqs. 16 and 17 in that paper, in that it yields both the 
redundant currents and the unknown voltage rises of current generators. 

We shall now show that Eqs. 3 are equivalent to the well-known 
Lagrange equations so widely used in solving dynamic mechanical 
problems. The latter equations can be typified as follows (1, p. 493): 


(8) 


where 7 is the kinetic energy of the system expressed as a function of 
certain selected parameters p2, Pm, called generalized coordinates, 
and their time derivatives (considered for purpose of differentiation as 
independent parameters). W is the net work done in increasing the 
kinetic energy of the system, or the difference between the work done by 
applied forces and that dissipated within the system or stored in it as 
elastic or potential energy or in any form other than kinetic energy. 
We shall postulate that inductive energy in the electrical network is 
analogous to kinetic energy in the mechanical system, thus: 


T = (9) 


so that inductance can be considered the analog of mass and current 
the analog of velocity, from which it can be shown that voltage and 
force must also be analogous. This corresponds to a well-known and 
widely used analogy between electrical and mechanical systems. 

By Eq. 1 the generalized coordinates in the electrical case are the 
prescribed and redundant charges q:, g2, ---@m- Hence the left side of 


(10) 


The right-hand side of Eq. 8 is called a generalized force, and is 
conventionally expressed as the increment AW of work done (other 
than kinetic energy) due to an infinitesimal increment Ap, divided by 
Ap. In the typical electrical branch jk the increment of charge due to 
Ag., say, is b;,,4g, and the voltage rise through which it is forced, with 
inductance ignored, is — — (Transformer voltage 
rises may be ignored because, taken together, they are not associated 
with energy, as can be proved in the particular case of Fig. 3 by sum- 
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ming voltage-current products, taking into account Eqs. 5.) The 
changes of voltage due to Ag, are ignored ; such changes must take place 
if g, is actually varied, but their omission corresponds to the fact that, 
in the conventional derivation of the Lagrange equations, the increment 
of the generalized coordinate p, is considered ‘‘virtual’”’ (4) in that it is 
assumed not to affect the forces within the system. Hence 


aw 
where the summation extends over all branches as before. By equating 
the right-hand sides of Eqs. 10 and 11 and comparing the result with 


Eq. 3, taking account of the definition (2), we see that the Lagrangian 
equations (8) are correct for the electrical systems under discussion. 


= [on — Riedie — (11) 


NODAL ANALYSIS 


Figure 4 shows two typical nodes in an electrical network, with 
constant values of capacitance, resistance and inductance arranged as 


N 
L jk 
44 
fi) 
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indicated. A prescribed current of instantaneous value i ;, is maintained 
by the indicated generator. The time-integral of voltage at the nodes 
(measured above some convenient datum) will be indicated as f; and f;, 
respectively, where f will for convenience be termed the flux (not to be 
confused with ordinary magnetic flux, although it has the same dimen- 
sions). As before, the parameters will be measured over and above the 
values at ¢ = 0. 

No matter how complicated the network, it is always possible to 
designate nodes so that each circuit element lies directly between two 


nodes as in Fig. 4. 
The total current flowing in the direction j to k can be shown to equal 


Cie (Fi— Fe) + (Fi — fe) / Rint (Fi — fie) /Lietin= Val tin, 


where Y,, is an appropriate operator. Since the total current leaving 
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node j must be zero, we may write 


CV fs — fe) + tin] = 0, 


where & is summed over all nodes other than j. In this expression 7 ;, 
may also be used to denote unknown current flowing through generators 
of prescribed flux rise, if any. 

If transformers exist, they can be represented by a combination of 
ideal transformers and appropriate circuit elements as before. Assume 
that Fig. 5, which is essentially the same as Fig. 3 with designations 
modified for our present purpose, represents the only transformer within 
a network. The fluxes at nodes 2 and 4 have been written in terms of 
those at 3 and 5 and the flux difference between 1 and h by consideration 
of the turns ratios. Applying the above equation to nodes h, 2 and 4 


in turn, we get: 


— fr) tin] = 0 
— fe) + tx] + = 0 
+ tar | + t' 45 0, 


where the summation terms extend over the entire network except for 
the transformer windings, and 7’ denotes transformer currents. If we 
multiply these equations by Nj:, Nes and N4s, respectively, and add 
them together, the sum of the transformer current-turns products equals 
zero by Eq. 4, and the remainder can be put in the following form : 


| [Vaal — fa) + + Se — fa) + 
+ Fi) + 0, (12) 


Al 


where the coefficient of each bracketed quantity is the coefficient of f, at 
the respective nodes in Fig. 5. 

We may generalize from this specific case as follows. Taking trans- 
former turns ratios into account, express the flux at any node 7 thus: 


fi = + + + (13) 


where the /’s on the right contain all prescribed flux differences and as 
many additional unknown node fluxes (to be termed redundant, like 
the redundant charges discussed previously) as are necessary to define 
each node flux. When no transformers exist the 6’s can equal only 
zero or plus or minus 1, whereas in the presence of transformers they 
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may have any constant values depending on the turns ratios. In any 
case the equations corresponding to Eq. 12 may be written as follows: 


fi — fr) + ta] = 0 s = 1,2,---n, (14) 


where transformer currents may be ignored ; the summation with respect 
to k is repeated for each node j in turn, and s can denote either a re- 
dundant node flux or a prescribed flux difference. The resulting equa- 
tions are just sufficient in number to yield the redundant fluxes and 
also the unknown currents associated with prescribed flux differences, 

To show that the redundants can be evaluated separately from 
the unknown currents of flux generators, denote each end of such a 
generator as a node, so that if the flux at one of these nodes is the re- 
dundant f, the flux at the other end is f, + f, or f, — f», where f, is 
the prescribed flux rise. In applying the above equation with s = r, 
the appropriate value of } will be unity at each node, and the current 
through the flux generator will be counted in a different direction at the 
two nodes, so that it disappears in the summation. 

The fact that transformer currents need not be considered explicitly 
may give considerable advantage over conventional methods of nodal 
analysis when transformers are present. If desired, these currents may 
be obtained by current-continuity considerations after the fluxes have 
been determined. 

In the steady-state case Eqs. 13 and 14 reduce to 


5 =1,2,---m}’ 
ik 


where Y is complex admittance. The second of these corresponds to 
Eqs. 19 and 20 in reference 3. 
To show that the Lagrange equations (8) apply to this case, we 
postulate that 
i sb 


where each branch is to be counted only once in the summation. Here 
capacitance is considered the analog of mass and voltage and analog of 
velocity, from which it can be shown that current and force must also 
be analogous. As in the previous case this corresponds to a familiar 
analogy between electrical and mechanical systems. 

By Eq. 13 the generalized coordinates in the electrical case are the 
prescribed and redundant fluxes. Hence 


— = bi Cul fr). (17) 
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To develop the generalized force corresponding to the right-hand 
side of Eq. 8, we express AW in terms of the incremental flux at each 
node due to the infinitesimal increment Af,, namely, 6;,4f,, and the total 
current approaching that node, where the incremental changes in cur- 
rents are ignored to correspond as before with the procedure used in the 
conventional derivation of the Lagrange equations. Hence, 


ow 
of. 
where transformer currents are ignored because, taken together, they 
are not associated with energy. By equating the right-hand sides of 
Eqs. 17 and 18 and comparing the result with Eq. 14, taking account 


of the symbolism of Y, we see that the Lagrange equations (8) are valid 
for the electrical network under discussion. 


ILLUSTRATIVE EXAMPLES 


bul (Si — fe)/Rie + (fi — t+ tn), (18) 


In Fig. 6 the voltage rise v,2 (where g designates ground) and the 
current i;, are prescribed, and the transformers’are ideal and their rela- 


Fic. 6. 


tive numbers of turns are indicated by expressing their voltage rises 
in terms of e, and é, respectively, the reference direction for the e’s 
being indicated with the aid of a plus sign placed at one end of each 
winding. The passive circuit elements are as indicated. It is desired 
to evaluate all branch currents, and all node voltages measured above 
ground. For simplicity all currents and voltages are considered zero 
= 0. 


Mesh Analysis 


The designations of flux should be ignored for the present. By 
applying Kirchhoff's first law establish that all branch charges can be 
written in terms of the prescribed charge (which is merely the time- 
integral of the prescribed current) and one redundant charge. Arbi- 
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trarily selecting the latter as g2:, and taking due account of the effect of 
the transformer turns ratios on the charges, we write all branch charges 
in terms of g74 and ge; as shown in Fig. 6. 

Next applying Eq. 3 with respect to g2; and gz, in turn, and ignoring 
transformer voltage rises, we have: 


1 on 2(- 2LG21) = 0 (19) 
— 1(— 042) + 1( — 074) — 2(— 2gr6/C) = OF 


As expected, the first of these suffices to determine the single re- 
dundant gs:, and reduces to: 


4L Goi + Roar = 


so that gz: may be evaluated if v,2 is specified and the resulting differen- 
tial equation can be integrated. The branch currents are then easily 
written by Kirchhoff’s first law in terms of 7;, and v,2, with results as 
shown in Fig. 7, where the symbolic operator D = d/dt is used in con- 
ventional algebraic manner. 

The second of Eqs. 19 yields the value of v;., which in conjunction 
with the voltage rises across the inductance, capacitance and resistances, 
easily calculable now that all charges are known, permits the voltage 
at all ungrounded nodes to be written, with results as shown in Fig. 7. 
As a check on the work the reader may evaluate the transformer voltage 
riseS V32, Vs36 and V5, Viz and see whether they obey the relationships 
assumed in Fig. 6. 


Nodal Analysis 


By consideration of the transformer turns ratios we may express 
all node fluxes in terms of the prescribed flux f,2: (which is merely the 
time-integral of the prescribed voltage) and the redundant node fluxes, 
namely, f:, fs and f;, as shown in Fig. 6, where the charge designations 
should be ignored for the present purpose. Applying Eq. 14 with 
respect to the three redundants in turn, and ignoring transformer 


currents, we have 
-1(fi/R) + 3(fi + fox)/2L = 0 
— fr)/R +4C(fo2 + fs)/2 = 0 
1(f: — frs)/R+in = 0. 
As expected, these equations do not contain the unknown current 
to: through the flux generator. The redundants may be expressed as 
follows: 


4f,/R fi/L fa2/L 
fi/R = fi/R — tn 


= 
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so that they can be evaluated in terms of the prescribed current and 
flux if the differential equations can be integrated. All node voltages 
are then easily written by consideration of the ratios of transformer 
voltage rises, with results as shown in Fig. 7. 

To evaluate i,2, apply Eq. 14 with respect to the prescribed flux 
difference, obtaining : 


1(— +3(fr + for)/2L + + fa)/2 = 0 


which yields the value shown in Fig. 7. The currents through the 
inductance, capacitance, and the resistance in branch 41 are found by 
considering the now-known voltage differences across these elements, 
and the transformer currents are determined by Kirchhoff’s first law. 
The entire solution checks that previously found by the mesh analysis. 


CHOICE BETWEEN MESH AND NODAL ANALYSES 


In complicated networks the difficulty of solution is determined 
largely by the labor of solving the differential equations defining the 
redundants, since the latter must be evaluated simultaneously. In turn 
this labor is dependent to a large extent on the order of the matrix of 
these equations, that is, on the number of redundants. Hence the 
choice between the mesh and the nodal analyses in any particular case 
should usually be in favor of that method which yields the fewer 
redundants. 

In the above example there was one redundant in the mesh analysis 
and three in the nodal analysis, so that the latter was more difficult to 
apply. The number of redundants in the nodal analysis could have 
been reduced to two, by considering the combination of the current 
generator and the resistance in branch 37 as a single new current gen- 
erator, with an over-all voltage rise v3, equal to the sum of the voltage 
rises in branches 37 and 74; this would have eliminated node 7 during 
the solution for the redundants. The voltage at 7 could easily have 
been found by separate calculation after v3, was determined. 

A corresponding procedure could have been used in the mesh 
analysis if the voltage generator were shunted by a passive circuit 
element, by considering the combination of the generator and the 
additional shunt element as a single new voltage generator, with a single 
over-all charge equal to — 
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ON A “BEST” FIVE-POINT DIFFERENCE ANALOGUE 
OF LAPLACE’S EQUATION 


BY 
DONALD GREENSPAN! 


1, INTRODUCTION 


A popular technique for the numerical solution of a Dirichlet prob- 
lem involves the replacement of the Laplace equation 


(1) 


Use + Uy = 0 


by the five-point difference analogue 
— 4uy + + Ue + U3 + = O, 


(2) 


where uo = u(x, = u(x +h, y), Us = u(x, y +h), us = u(x —h, y), 
us, = u(x, y — h), and h is the positive grid constant. (See, for ex- 
ample, reference (1),? chapters VIII and X, for a complete discussion.) 
Also, it is known that Eq. 2 is a “‘best”’ five-point analogue of Eq. 1, 
subject to certain assumptions and definitions. (See, for example, (2), 
or the works of Collatz.) 

The object here is to generalize the known results and to demonstrate 
a ‘“‘best”’ five-point analogue of Eq. 1 for which the grid size h, in the x 
direction, need not necessarily be the same as the grid size d, in the y 
direction. 


2. DERIVATION OF THE DIFFERENCE ANALOGUE AND REMARKS ABOUT IT 


Let the grid sizes in the x and y directions be / and d, respectively 
(see Fig. 1). The following notation will be used throughout: 


uo = u(x, y), = u(x +h,y), us = u(x,y +), (3) 
u;= u(x —h,y), us = u(x, y — 


ox™dy" 


= Una; = Uma (4) 


Since 4 and d are assumed to be positive constants, let 


h = pd 


(5) 


where ? is a positive constant. 


1 Mathematics Department, Purdue University, Lafayette, Ind. 
2 The boldface numbers in parentheses refer to the references appended to this paper. 
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Since u is harmonic, then tn42,. + Um,n42 = 0; m, m non-negative 
integers. Hence: 


wu 


In terms of the A’s, (6) may be rewritten, by use of (5), as follows: 


(2) 2Aiat Az2=0 
(j) 10A 0,5 + = 0 
(Rk) 15A6¢,0 + Aas = 0 
(2) 10A5,1 + 3A3,3 0 
(m) =0 
(nm) 10A;,5 + 3A33 = 0 
(0) 15Ao6e+ = 9, 


6A40 + 
Asi 
6Ao,1 + 
10As0 
2A41 + 


40 [J. F. 1. 

| 
(xth,4) x 
3 0} (<4) 
4 Cx, 
(a) Ue, o + Uo,2 = 0 (2) U3,2 + U1,4 
(b) U3,0 + 412 0 U2, 3 + Ko, 5 
(c) Ue,1 + = 0 (R) U6,0 + U4,2 
(d) Uso + = 0 (1) + (6) 
(ec) (m) Use + 
(f) to. = 0 (mn) Uss+ 
(g) UsotuUs2=0 (0) Us4+ 
(h) t+ =0 
(a) Asoo + = 

3A30+Ai2=90 

(c) 3Ao3 

(d) 7 
(e) 
(f) =0 
(g) 
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or, equivalently, 


(a) Aso Ao,2 (4) = 2Ai4 


(6) Aso = — (1/3)Ai2 (j) — 

(c) Az1 = — 3Ao,s (k) Aso = — Ane 

(d) Ags As1 =Ais (8) 
(e-) Asi = —Ais (m) Ags = 15A06 

(f) = — Ass = — (10/3)Ai5 

(g) Aso = (1/5)A14 (0) Aas = — 15Ao.. 


(h) 


Now let 


Agi 5Ao,5 


L(u) = (9) 


For future reference, recall by Taylor series about (x, y), and by use of 
(5), that: 


uy = Ao,o 


Ao,o + Ai, opd + A>, op*d? + A;,op'd® + A, op*d* 
+ As, op'd® + Ag op*d® + O(a’) 


= + Ao id + Ao,.d* + Ao, sd' + Ao, 
+ Ao, + Ao, d* + O(d’) (10) 


us = — Aropd + Ax op’d? — Az, op*d* + Ax 
— As op'd® + Ag op'd* + 0(d") 


= Ao,o Ao, id + Ao? — Ay Ao, 
— Ao, + Ao + O(a’). 


Substitution of (10) into (9) and recombination yields: 


L(u) = + a1 + a2 + a3 + a4) + PdAi,o(a1 — as) 
+ dAo,1(@2 — as) + p*d*Ae o(a1 + Gs) + d*Ao, + a4) 
+ — a3) + d*Ao — as) + o(ai1 + 
+ d'A + a4) + 5, a3) + as) 

+ p%d*A + a3) + d*Ao,6(@2 + a4) + O(d"). (11) 


0 
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By use of (8), (11) becomes 


L(u) = Ao,o(@o + a: + G2 + as + G4) 
d[A 1 op (a1 a3) +f ~ a.) | 
+ (a2 + a4) + as) 
+ d*[Ao,s(@2 — a4) — (1/3)p*A1,2(@1 — as) 
+ + as) + (a2 + a.) | 
+ d[Ao,s(a2 — as) + (1/5)p°A1,4(a1 — @s)] 
d°Ao, (a2 as) p*(ai + a) | + O(d"). (12) 


In order to make (12) of as high an order in d as possible, consider 
the following set of equations 


a) +a, + 42+ 4; + a = = Od) 
a; — = = 0(d*) 
— dy = = O(d*) 
G2 + — + G3) = & = O(d*). 


(12’) 


The solution of this system in terms of dp is: 


—2(1+") —2(1+") 
(13) 
_ testprertes _ — +estesp?— ea 


Hence, (9) becomes, finally 


— — eof? — + + ao — — — — + a 


do — €1 + + peo + &% — + + — 


U4. 


The choice of €2, is governed by an attempt to make the 
right hand side of (11) of as high an order ind as possible. If one chooses 
ay = — 4, €1 = 2 = €; = & = O, then a; = a; = 2/(1 + p*), a2 = a 
= (2p*)/(1 + p*). Hence (14) becomes 


(U2 + Us) (15) 


L(u) = — + + + 


while (12) becomes 


L(u) = 0(d*). (16) 


; 
‘ 
: 
(14) 
14 
: 
t 
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From (15) and (16) then 


+ (uy Us) + —— (us U4) 0(d‘) = = 0. (17) 


Discarding the terms 0(d*), one readily deduces from (17) the desired 

difference equation — of (1), that is, 

(u2 + U4) = 0. (18) 


— 4u) + —— + U3) + 


Note that for the special case p = 1, that is, when h = d, (18) re- 
duces to (2). It should also be noted that (18) has been derived else- 
where by other methods, but that the techniques used here will prove 
to be of great value in what follows. 


3. ESTABLISHMENT OF EQ. 18 AS A “BEST” ANALOGUE,OF LAPLACE’S EQUATION 


In the discussion of Section 2, Eq. 17 is the key equation. However, 
it must, in general, be studied carefully. Suppose, for example, one 
had set a) = — 4d? and e; = é: = e; = ¢s = 0 in (14). Then (17) 
would have taken the form: 


+ —— (t1 Us) Us) 0(d*) = 0. (17’) 


1 2 oe 
The process analogous to that described in going from (17) to (18) 
would yield in the case of (17’) 


+ (uy 3) +. (tue Us) = 0. (18") 


It appears that (18’) results from (17’) by discarding all terms of 
order at least 6 ind. But, of course, (18’) can be divided by d?, so that 
this latter statement is correct, but misleading. Such trouble is avoided 
as follows. A difference analogue of (1) will be called allowable if and 
only if 


4 
(a) it is of the form: ¥ ama; = 0, each u; defined by (3), and 
0 


(6) lim a; ¥ 0, for at least one value of i = 0, 1, 2, 3, 4. 
Suppose then that in general one obtains in place of (17): 


> au; + 0(d") =0 (17”’) 


0 
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and in place of (18): 


4 
au; = 0, (18”’) 
0 


and suppose that (18’’) is an allowable difference equation, then (18’’) 
is called an approximation of order (” — 1). 

It follows immediately that (18) is an allowable third order ap- 
proximation. 

Now it will be shown that (18) is a “best” difference equation 
analogue of Laplace’s equation in that both the following conditions 
are satisfied : 


(a) there does not exist an allowable fourth order analogue, and 
(6) every other third order analogue differs from (18) by only a 
multiplicative constant. 


Theorems 1 and 2, below, establish these results. 

Theorem 1. There does not exist an allowable fourth order difference 
equation analogue for Laplace’s equation. 

Proof. Suppose such an allowable fourth order approximation 
exists. It follows from (a), (8), (c), (d), (e), (f), of Eqs. 8, that Ao >, 
A;,2, Ao,s, A1,3, Ao,4 may be viewed as independent parameters for this 
sub-set of equations, while A», A2,2, As,1, may be con- 
sidered dependent variables. Moreover, since Ao,o, Ao,1, A1,0 are inde- 
pendent of d, it follows that: Ao,o, Ao,1, A1,0, Ao,2 A1,2, Ao,s, A1,3, Ao,4 
are all independent parameters. Hence, if the right sides of (11), and 
hence (12), are to be of order at least 5 in d, then it must be that 


ao +4, + +43 + = = Od) 

a; — a; = = 

G2 — = = O(d*) 

(a2 + a4) — + a) e, = 0(d*) 
+ as) + a2 + ay = es = O(@') = Od). 


The solution of this sytem is 


p(itp) 1+p*’ 2 °2(1+p*)p” 


Hence, lim a; = 0, for each ¢ = 0, 1, 2, 3, 4, by definition of the e;, 
d—0 


i = 0, 1, 2, 3,4. But this is a contradiction to the assumption that an 
allowable fourth order analogue exists. Hence the theorem follows. 


a2 


U. 
af 
2 
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Theorem 2. There does not exist an essentially different, allowable 
third order difference analogue of Laplace’s equation other than Eq. 18, 
that is, every other allowable third order analogue differs from (18) only 


by a multiplicative constant. 

Proof. The discussion of Section 2, through the development of 
Eq. 14, is that which is sufficient to establish the general third order 
analogue. Hence consider Eq. 14, rewritten as follows 


X[ei(— — Uy — Usp”) +e2(— Ui — Usp? +u3+usp’) 
+ +64 (U1 Uy) |. (19) 
By definition of ¢:, 2, €3, ¢, in (12’), and by use of (10), it follows that 


usp”) = 0(d*) 

€2(— U1 = =0(d*)0(d) (20) 
U2+U3— U4) =0(d?)0(d*) =0(d*). 


By use of (20), (19) may be rewritten: 


In the same way that (18) was established, it follows that the gen- 
eral, allowable third order analogue is given by 


+ Us _ + us) 


where a» is an arbitrary, non-zero constant. Hence, from Eq. 22, the 
theorem follows. 
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Model Magnified Billion Times.— 
A model of a molecule magnified one 
billion times was fabricated by an Air 
Research and Development Command 
center for exhibition at the 1958 
Brussels World’s Fair. It is a model 
of ferrocene, an American-discov- 
ered chemical compound. ARDC’s 
Wright Air Development Center is 
studying derivatives of ferrocene as 
fluids and lubricants for space vehicles. 

To convert ferrocene into high- 
boiling liquids, WADC scientists dis- 
covered that the compound, which 
has unusually good resistance to high- 
temperatures, can be modified to 
provide lubricant characteristics and 
a wide liquid range. In early tests 
of ferrocene compounds as potential 
lubricants, one derivative of the com- 
pound retained its fluidity throughout 
a 1000 degree temperature range from 
minus 112 degrees to more than 900 
degrees Fahrenheit. 

Ferrocene was discovered in 1951 
by two chemists, T. J. Kealy and P. 
L. Pauson of Duquesne University, 
who came upon it while experimenting 
with another compound. It was not 
until several years later when the 
compound’s molecular structure was 
determined that more extensive re- 
search was begun. 

Ferrocene, which is an organomet- 
allic compound, gets its unusual ther- 
mal stability from its molecular 
arrangement. A “sandwich” com- 
pound, considered rare by chemists, 
ferrocene has proved to be the most 
stable among the “metal-sandwich”’ 
structural types. This unusual struc- 
ture also makes possible the first 
known chemical compound composed 
of carbon, hydrogen and iron. 

The name ferrocene is derived from 
ferro, for iron, and cene, for benzene 
which it resembles in chemical 
reactivity. 


CurrRENT Topics 


{J. F. 


Another outstanding property of 
ferrocene is its anti-soot effect. In 
tests using ferrocene as an additive to 
fuel oil, early research workers found 
no soot or carbon formed in the fire 
chamber even after prolonged periods 
of burning. 

A ferrocene additive has been in- 
vestigated for use in jet-engine fuel 
and oil to eliminate carbon deposits. 
A compound related to ferrocene is 
being considered as an anti-knock 
additive for automobile engines. 

WADC’s Dr. Harold Rosenberg 
planned the ferrocene exhibit as a 
member of the U. S. Molecule (Chem- 
istry) Committee for U. S. participa- 
tion in the Fair. He is a senior proj- 
ect engineer in the Organic Materials 
Branch of the Materials Laboratory 
at WADC. 


New Device Simplifies Transistor 
Mounting.—Transipad, a new device 
which simplifies installation and per- 
mits a more reliable mounting for 
transistors on printed circuit boards, 
has been announced by Milton Ross 
Metals Company, Southampton, Pa. 

Designed especially for JETEC 30 
transistors, the Transipad is a glass 
filled Diallyl Phthalate wafer with 
three holes and three hemispherical 
feet. The transistor leads pass 
through the holes to the circuit board 
solder connections, and the transistor 
rests directly on the Transipad. This 
eliminates the necessity of using the 
leads as supports and permits a lower, 
more stable mounting with positive 
insulation between the transistor case 
and the printed circuit conductors. 

The hemispherical feet raise the 
Transipad off the board, allowing flow 
clearance for the solder fillets and 
permitting air to circulate between the 
circuit board and the Transipad. 
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THE COSMOLOGICAL PRINCIPLE AND THE 
COSMOLOGICAL CONSTANT 


BY 
PARRY MOON! AND DOMINA EBERLE SPENCER? 


1. INTRODUCTION 


Cosmological theories, from Einstein’s work in 1917 to the recent 
models of Bondi and Hoyle, place strong reliance on the so-called 
“cosmological principle.” This is, of course, not a principle of nature; 
‘it is a mere hypothesis. The idea, however, is a reasonable one: since 
we have almost no experimental facts about the cosmos, let us make 
the simplest general assumption that we can imagine, namely, that the 
universe appears the same to all observers. 

Obviously, this hypothesis is based on little or no direct evidence : 
indeed, most astronomical observations seem at first sight to deny its 
validity. The familiar constellations are present only for an observer 
in the vicinity of Sol. The configuration of the nearby galaxies would 
be replaced by another configuration, if the observer moved to a different 
galaxy; and even the clusters of galaxies at the limit of telescopic 
observation would have a different appearance to a distant observer. 
If one accepts Mach’s principle, this variation in the surroundings may 
be expected to result in some variation in the laws of nature as measured 
in different parts of the cosmos. At least, we have no proof that it does 
not. Thus the cosmological hypothesis, if true at all, has validity 
only in a statistical sense; and the exact statement of this statistical 
assumption is not easy. 

Various definitions have been formulated in attempts to obtain a 
rigorous statement. Einstein said (1),* “‘Alle Stellen des Universums 
sind gleichwertig.”” Milne obtained greater precision by introducing the 
idea of a substratum and an associated set of equivalent observers. 
Two observers are said to be equivalent (2) ‘‘when the totality of ob- 
servations A can make on B can be described by A in the same way as 
the totality of observations B can make on A can be described by B.” 
Bondi (3) says, ‘““The cosmological principle may now be stated as: 
Every member of a three-parameter set of fundamental observers ob- 
tains the same result of corresponding observations of the universe as 
every other member.” 

One may introduce the additional assumption that the cosmos is in 

1 Department of Electrical Engineering, Massachusetts Institute of Technology, Cam- 
bridge, Mass. 

2 Department of Mathematics, University of Connecticut, Storrs, Conn. 

3 The boldface numbers in parentheses refer to the references appended to this paper. 
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statistical equilibrium, so that time as well as position has no effect. 
Bondi calls this the ‘perfect cosmological principle.”” With either the 
ordinary or the “perfect’’ hypothesis, the introduction of an idealized 
substratum lessens vagueness and allows the building up of a deductive 
science, as has been done by Milne (2). The weakness of this approach 
seems to lie in its divorce from reality. Everything depends on hypo- 
thetical measurements of idealized, non-existent observers; so that one 
may have an uneasy feeling that the whole deductive system has little 
relation to the real world. 

A third criticism of the cosmological hypothesis is that it has not 
justified the hopes of its originators. Many models of the universe 
have been devised, but none is satisfactory. To summarize, the cosmo- 
logical hypothesis may be criticized because 


1. It has no direct experimental justification. 

2. The difficulty of obtaining a rigorous definition leads to the 
suspicion that the cosmological hypothesis has little relation to the 
actual universe. 

3. As a tool for the creation of a unique cosmological model, it has 
been unsuccessful. In fact, as long as we maintain the fundamental 
requirement that the cosmological hypothesis be satisfied, we are 
raising an artificial barrier against further progress. 


The cosmological hypothesis was introduced as a temporary ex- 
pedient in view of the dearth of experimental data. But a considerable 
amount of information has accumulated during the past few years, so 
that we are now in a much better position than Einstein was in 1917. 
The time is ripe for an inquiry into the possibility of abandoning the 
cosmological hypothesis and basing cosmology on experimental facts. 


2. THE NEWTONIAN APPROACH 


The subject of cosmology is usually associated with the names of 
Einstein, de Sitter, Eddington, and Lemaitre. Their models are based 
on general relativity and therefore involve a form of the cosmological 
hypothesis. The original Einstein model was unstable as well as non- 
expanding. The expanding model of de Sitter contained no matter and 
has therefore only the remotest relation to the actual universe. Further 
modifications by Eddington, Lemaitre, Hoyle, and others have failed 
to produce a unique model that satisfies the requirements. 

Certainly, a satisfactory cosmology must not contradict established 
principles of terrestrial physics, though it may indicate that these 
terrestrial relations are only special cases of more general laws. Having 
this idea in mind, we experience a distinct shock when we discover that 
the expanding relativistic models specify infinite velocities. Experience 
with particle accelerators shows that the velocity of electrons with 
respect to the accelerator can never exceed the velocity of light. If the 
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velocity of elementary particles has an upper limit it is likely that the 
velocity of galaxies is likewise limited. 

For this and other reasons, we tentatively reject cosmologies based 
on general relativity. The physicist and the engineer employ New- 
tonian mechanics in dealing with all ordinary macroscopic phenomena ; 
and it is only reasonable to employ these familiar relations as the basis 
for a large-scale, cosmological formulation, as has been done by Milne 
and McCrea (4). In taking this step from the moderately large to the 
very large, we may expect that some modification will be necessary in 
Newton’s equations; but the more general equations should reduce to 
Newton’s for moderate distances. Indeed, one of the principal fruits 
of cosmological study may well be the additional light it can throw 
on the equations of every-day physics. 

We propose to modify Newton’s gravitational equation by the 
inclusion of a factor : 


(1) 


= #(£), 


where 


#2: = force on a concentrated mass m:, caused by mass mi, 
a: = unit vector pointing from m, to mz, 
& = distance between m, and mz, 


#(£) = a gravitational function having the value —1 for small distances 
and becoming positive for galactic distances. 


Evidently some such modification must be made in cosmology to account 
for the recession of the galaxies. 

Equation 1 applies to two mass-points. For a continuous distribu- 
tion of matter of density p, the force on m, becomes 


= mG f (2) 


For the special case of a symmetric, spherical universe with p = p(t) 
and # = const, 


(2a) 


where r is the distance from m, to the center O of the cosmos. Equation 
2a, obtained by direct integration, agrees with the usual theory of 
gravitational potential. The shell of matter outside r has no effect on 
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m,, and the force is as if the total mass of a sphere of radius r were 
concentrated at O. This simple result, however, is obtained only if the 
density is uniform and @ = const. 

In the Newtonian formulation of cosmology, we use the equation of 
motion and the equation for the conservation of matter: 


(3) 


(4) 


div (pv) 


It is customary to employ a spherically symmetric model of finite radius 
in euclidean 3-space. Velocity is radial and p is a function of time only. 
Then the foregoing equations become 


= Go) f d&, (3a) 


Solution of Eqs. 3a and 4a gives the density p(t) and the velocity v(r, ¢). 
It is customary to assume that the universe is expanding, with 
velocity directly proportional to radius: 


v = F(t)-r. (5) 


Substitution into Eqs. 3a and 4a yields 


= 


If 6 = const, elimination of p results in the differential equation, 


dF 

Note that Eq. 6 is independent of the value of @ = const. Thus if 
one assumes linearity, Eq. 5, the velocity distribution is given directly 
by a solution of Eq. 6, irrespective of ®. 


3. EXPANDING MODEL, MILNE I 


As a simple example, consider a modification of the model de- 
veloped by Milne (4) in 1934. The universe consists initially of tightly 


dv a:p®(é) 
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packed material in a sphere of radius 7); and this sphere explodes at 
t = 0, sending out fragments having all velocities from 0 to c.4 The 
center of the universe is taken as the origin O of a spherical coordinate 


system. 
Assume that the gravitational force is expressed by the usual 


Newtonian formula. Then @ = — 1 and the right side of Eq. 3a 
becomes 


Gp (é)r, 


in accordance with Eq. 2a. A solution of the problem, obtained from 
Egs. 6, 3a, and 4a, is 


(7) 


where a is an arbitrary constant. 


At any instant ¢, the total mass M of the universe is contained in a 
sphere of radius 7,(#). Thus 


M = 


E (t + (9) 


The maximum velocity, at a given ¢, occurs at r = 7, (t); and from 


E uc)" (t + a)-™, (10) 


The initial condition is now introduced. Att = 0, Umax = c¢. Thus, 
from Eqs. 10 and 9, 


4MG 


3c8 ’ 


(11) 


2MG 


(12) 


= 


4 Milne assumed that velocities ranged from 0 to « and that ro = 0. 


; 
= 
3t+a’ 
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: Eqs. 7 and 9, 
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Each galaxy is undergoing continuous deceleration, and 


9 (¢ + a)* 


The behavior of the model is indicated in Figs. 1 to 3. 


c 


t 


Fic. 1. Radius of the cosmos, according to Fic. 2. Maximum velocity of galaxies 
the expanding model, Milne I. (Milne 1). 


li 


Fic. 3. The velocity distribution is a F1cG.4. Radius of the cosmos according to the 
linear function of radius, but the slope of the unaccelerated expanding model (Milne II). 
line decreases with time (Milne I). 


4. FINITE OR INFINITE? 


Note that the simple model of Section 3 is obtained without any 
modification in the ordinary Newtonian equations. Also, it satisfies 
Olbers’ principle, galactic velocities are directly proportional to dis- 
tances, and maximum velocity does not exceed c. 

But, strangely enough, the model fails to show the one thing that 
Milne was trying to prove, namely, that the models of general relativity 
can be obtained exactly by Newtonian methods. He says (4), ‘‘In this 
paper I show how the same locally observable results can be obtained 
from elementary Newtonian theory as are given by Einstein and de 
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Sitter’s well-known relativistic model. . . .”’ Also, “‘. . . these [equa- 
tions] are identical with the relativistic equations of an expanding 
universe of zero curvature with pressure = 0 and cosmical constant 
= 0, as given by Einstein and de Sitter.’’ This conclusion has been 
generally accepted. Actually, Milne I is not the Einstein-de Sitter 
universe, nor does it satisfy the cosmological hypothesis. 

A basic characteristic of relativistic cosmologies is that every point 
is equivalent to every other point. But the Milne models have a 
preferred point, the origin O from which the expansion takes place, and 
this point is uniquely distinct from all other points. One can easily 
prove (5) that, for any system satisfying the linear relation, Eq. 5, each 
observer will see the other galaxies moving away from him on all sides 
and will observe a statistical isotropy. But this is true only if the 
limit of his telescope is less than the apparent distance to the boundary 
of the universe. And the Milne universes do have a sharply defined 
boundary, beyond which there is nothing. An observer near the edge 
of the universe will certainly not obtain the same picture as an ob- 
server at O. In other words, the cosmological hypothesis is not valid, 
except in a very restricted sense. 

A similar criticism of the Milne derivation was proposed by Layzer 
(6). Inreply, McCrea (5) claimed that though the derivation was made 
for a finite sphere, the radius 7, can be allowed to expand to infinity; in 
which case, the cosmological hypothesis holds. But if r; = ©, sym- 
metry and homogeneity require that the gravitational force be zero 
at any point. Then the right side of Eq. 3a becomes zero, and Eqs. 7 
and 8 are no longer solutions. Also, if matter is conserved and if the 
total mass M (of the original universe of radius 7») is finite, then the 
density p of the infinite universe must be zero. The solution given in 
Section 3 applies only to a finite, spherical, expanding universe; and 
the cosmological hypothesis does not hold, except as an approximation 
over a limited region. 

This conclusion should not be construed as a reflection on the im- 
portance of the work of Milne and McCrea (4). In fact, their discovery 
of the Newtonian technique in cosmology seems of great value, particu- 
larly in view of the unfruitfulness of the relativistic approach. 


5. EXPANDING MODEL, MILNE I 


Milne I was obtained without any alteration in the Newtonian 
equations. Now consider a more general treatment, where the gravita- 
tional force is represented by Krp(t) on the right side of Eq. 3a. If Eq. 5 
holds, the differential Eq. 6 still applies, irrespective of K. 

As a special case, we now take the unaccelerated model used by Milne 
in nearly all his work (7). Here K = 0, so Eq. 3a becomes 


dv dv 


(30) 
or at 
« 
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with a solution 


r 
(14) 


Suppose that the universe starts expanding from an original radius fo 
at ¢ = 0, and that the maximum velocity isc. Then 


(15) 


the radius of the universe is 


(16) 


and the velocity is 


r 
(14a) 


The total mass is 


M = p(t) 


3M 
A4nc®(t + ro/c)* (17) 


Characteristics are shown in Figs. 4 and 5. Particularly simple results 
are obtained by setting ro = 0, as done by Milne (4). 


VA 


Fic. 5. The velocity is always a linear function of radius. The maximum velocity 
always occurs at rf =7; and is always c. (Milne IT). 


Various other values of K may be employed to obtain expanding 
universes with positive or with negative accelerations. But in all cases, 
the universe is finite and the cosmological hypothesis is not satisfied. 
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6. THE COSMOLOGICAL CONSTANT 


In this family of expanding models, Sections 3 and 5, the gravita- 
tional force is Krp(t) and K may be written 


(18) 


The quantity \ was originated by Einstein in 1917 and is called the 
cosmological constant. Its purpose is to introduce an ad hoc force of 
repulsion and thereby allow an expanding universe. In Milne I, the 
universe expands because of the initial velocity of its particles, despite 
the fact that \ = 0. In the other models (Section 5), however, \ > 0. 
In Milne II, K = 0, and thus \ must be adjusted so that for every point 
and for every instant, the two terms of Eq. 18 exactly cancel. 

The idea of the cosmological constant has always been a disturbing 
one. In fact, the concept seems at least as repugnant and as hard to 
visualize as the concept of an aether. We shall attempt to show, 
however, that an associating \ with a modification in Newtonian gravi- 
tation, a reasonable picture is obtained. 

If K = 0, for instance, the Newtonian gravitational force must be 
cancelled completely over the entire universe. But we know that 
Newtonian attraction occurs in the solar system and even within 
galactic clusters, though repulsion may occur between clusters. This 
inconsistency—the rigid requirement that K = 0 everywhere, con- 
trasted with the experimental fact that K is negative for moderate 
distances—is a result of the usual treatment of the cosmos as a con- 
tinuum or as a sort of fluid. If the discrete nature of the galaxies is 
taken into account, the inconsistency disappears. 

Our Eq. 1 allows galactic recession by introducing the function 
#(£), which must be negative at small distances and positive at large 

‘distances. Data are not yet sufficient to allow an exact specification of 

(£), though one may assume tentatively that 


= (19) 


where 6 is the distance at which ® = 0, and ®,, is the value of ® at 
r-—»«, This function is indicated in Fig. 6. 

If 6 be taken as2 X 10° lightyears, the effect of 6(£) on the behavior 
of planetary systems, binary stars, and even stellar motions within a 
galaxy, is negligible. Also, gravitational forces within the local group 
will still be attractive, tending to consolidate the cluster; but the forces 
between galactic clusters will be repulsive because of the greater dis- 


{J. F. 
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tances. For instance, if ) = 2 K 10° and #, = + 1, and we consider 
the gravitational effect between the local group and the nearby Virgo 
cluster (£ = 15 X 10*ly), Eq. 19 gives a repulsive force with a nu- 
merical value of 77 per cent of the Newtonian value. For more-distant 
galaxies, the repulsion will approach 100 per cent of the magnitude of 
the Newtonian force. 

If Eq. 1 be applied as if the cosmos were a fluid, no satisfactory model 
seems possible. But if we realize that galaxies are discrete, with large 


5 
2 
© 
2 = 
c 
| 
| 
10° io? io? 


Fic. 6. Proposed gravitational function of @(¢) for #, = + 1, giving an attractive 
force for — < 6 and a repulsive force for — > 6. 


spaces between clusters, a simple treatment is obtained. Because the 
distances involved are large compared with 6, @ = ®, and the integral 
of Eq. 3a reduces to a summation of terms of the form 


+a (1a) 


Consider a galaxy at point P, distant r from 0. The inverse-square 
relation of Eq. 1a allows the neglection of all galaxies beyond radius r, 
as in ordinary Newtonian potential theory. Thus the total force at P 


By 
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is outward, caused by all the matter within a sphere of radius 7, and 
Eq. 1a becomes 


dv ov 

The constant K, which was employed in Section 5, is expressible 
in terms of the arbitrary constant ©, : 


v + Gor®,,. (3c) 


Ge.. (20) 
Also, from Eq. 18, 


For Milne I, 


For Milne II, 


By choosing any other value of ®,, we obtain another model, but all 
of them are alike in not satisfying the cosmological hypothesis. We 
have purposely omitted models with continuous creation of matter, 
though they can also be handled by the Newtonian method. 


7. CONCLUSIONS 


Forty years of intensive work with general relativity have failed to 
produce a satisfactory cosmology. This fact raises a suspicion that 
the cosmological hypothesis is at fault—that a better approach to 
cosmology is a modification of the Newtonian method aided by all 
experimental facts now available. Conclusions may be summarized 
as follows: 


1. The cosmological hypothesis should be rejected. 

2. Newtonian models do not satisfy the cosmological hypothesis, 
despite frequent statements to the contrary in the literature. 

3. Newtonian models are not equivalent to relativistic models. 

4. For very large astronomical distances, Newton's gravitational 
equation needs revision. A proposed change in the equation allows 
the usual Newtonian attraction within galactic clusters but provides a 
possible repulsion between clusters to allow an expanding universe. 


#(£) = —1, 
4rG 
K=- = 
1 = 0. 
$(&) = 0, 
K = 0, 
4rGp 
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COMMITTEE ON SCIENCE AND THE ARTS 


(Abstract of Proceedings of Stated Meeting held Wednesday, June 11, 1958.) 
HALL OF THE COMMITTEE, 
PHILADELPHIA, JUNE 11, 1958. 


Mr. JoserH Gray JACKSON in the Chair. 


The following reports were presented for final action: 


No. 3214: The Atomic Clock. 

This report recommended the award of a Certificate of Merit to Harold Lyons, of Culver 
City, California, ‘In consideration of his work in pioneering the development of clocks of very 
great precision based on the natural periods of vibration of atoms and employing microwave 
techniques.” 

No. 3286: Rosenberg’s Cross-Field Dynamo. 

This report recommended the award of a Howard N. Potts Medal to Emanuel Rosenberg, 
of Bogota, Colombia, ‘‘For his conception of a fundamentally new principle for deriving elec- 
trical current from a rotating machine and for his solution of the many technical problems 
related to the development of the Cross-Field Generator.” 


No. 3313: Work of William N. Goodwin, Jr. 

This report recommended the award of a Howard N. Potts Medal to William Nelson 
Goodwin, Jr., of New Hope, Pennsylvania, ‘In consideration of his many contributions to the 
field of electrical measuring instruments of which the thermal ammeter, the vacuum tube 
tester, and the rectifier voltmeter are examples, and in particular for his work on photoelectric 
exposure meters and establishing a scale for film speeds.” 

No. 3323: Work of George S. Crampton on the Development of the “‘Borescope.”’ 

This report recommended the award of an Edward Longstreth Medal to George S. 
Crampton, of Philadelphia, Pennsylvania, “In consideration of his development of the ‘Bore- 
scope,’ a specialized periscope for the internal inspection of turbine shafts, rotors, guns, and 
various closed vessels, by which flaws may be detected thus providing protection of vital im- 
portance from the physical and economic hazards of operational failure.” 


No. 3343: Clamer Medal. 

This report recommended the award of the Francis J. Clamer Medal to Julian Miles 
Avery, of Framingham, Massachusetts, “In consideration of his discovery and application of 
the principles of high top-pressure operation of blast furnaces which have very materially 
increased their production and contributed to their smoother and more economical operation.” 


No. 3344: Brown Medal. 

This report recommended the award of the Frank P. Brown Medal to Charles Milton 
Spofford, of Brookline, Massachusetts, ‘For the engineering, aesthetic and educational ac- 
complishments of his life work, exemplified in his prompt recognition, espousal and clarification 
of valid structural theories ; his design of many large and beautiful bridges, and other prominent 
engineering works ; his contributions to engineering education through teaching and administra- 
tion of an engineering educational department ; and his authorship of widely used textbooks.” 


No. 3346: Henderson Medal. 

This report recommended the award of the George R. Henderson Medal to Arthur H. 
Morey, of Erie, Pennsylvania, “In consideration of his engineering contributions to, and his 
effective leadership in, the over-all design, development and service testing of the successful 
4800 h.p. General Electric Company Gas Turbine-Electric Locomotives which involved the 
close and skillful coordination of the specialized developments by the Turbine, and other com- 
ponent divisions of the General Electric Company.” 


D. S. FAHRNEY, 
Secretary to Committce 
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The Committee on Library desires to add to the collections any technical works that 
members may wish to contribute. Contributions will be gratefully acknowledged and placed 
in the Library. Duplicates received will be transferred to other libraries as gifts of the donor. 
583 items have been added the past month. 

Photostat service. Photostat prints cf any material in the collections can be supplied 
on request. 

The Library and reading room are open from 9 A.M. until 5 P.M. on Mondays, Tuesdays, 
Thursdays, and Fridays; from 2 P.M. until 10 P.M. on Wednesdays; and from 9 a.m. until Noon 
on Saturdays. 
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FRASER, RONALD GEORGE JuTA. Once Round the Sun. 1957. 
Moore, Patrick A. Guide to Mars. 1957. 
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NATIONAL Book LEaGuE. Science for All; an Annotated Reading List for the Non-Specialist. 
1958. 


BIOLOGICAL CHEMISTRY 


BREDEMANN, Gustav. Biochemie und Physiologie des Fluors und der Industriellen Fluor- 
Rauchschaden. Ed. 2. 1956. 
Moore, THoMAs. Vitamin A. 1957. 


BOTANY 
LeopoLp, ALpo Cart. Auxins and Plant Growth. 1955. 
CHEMISTRY AND CHEMICAL TECHNOLOGY 


Advances in Protein Chemistry. Vol. 12. 1957. 

AMERICAN CHEMICAL SociETY. Division oF PETROLEUM CHEMISTRY. Techniques of Cata- 
lyst Preparation. 1956. 

Arma, K. AND OTHERS, Yeasts. 1957. 

BERLow, EVELYN; BARTH, ROBERT H. AND SNow, JoHN E. The Pentaerythritols. 1958. 

BIKERMAN, JACOB JOSEPH. Surface Chemistry. Ed. 2. 1958. 

BurriEL Marti, FERNANDO AND RAM{fREz-MuNoz, J. Flame Photometry. 1957. 

HARNED, HERBERT SPENCER AND OWEN, BENTON B. The Physical Chemistry of Electrolytic 
SoLutions. Ed. 3. 1958. 

LANGHAMMER, GUNTER. Versuche zur Physikalischer Chemie. 1956. 

Puptett, P. A. R. Synthetic Detergents. 1957. 
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SmoTHERs, W. J. AND CHIANG, Yao. Differential Thermal Analysis. 1958. 

Watt, GEORGE WILLARD. Basic Concepts in Chemistry. 1958. 


DICTIONARIES 
Cassell’s German Dictionary, German-English, English-German, Rev. Ed. 1958. 
ELECTRICITY AND ELECTRICAL ENGINEERING 
BENSON, FRANK ATKINSON. Voltage Stabilized Supplies. 1957. 
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Programs for an Electronic Digital Computer. Ed. 2. 1957. 


ENGINEERING 


INTERNATIONAL CONFERENCE ON SOIL MECHANICS AND FOUNDATION ENGINEERING. 4TH, 
Lonpon, 1957. Proceedings. Vols. 1-2. 1957. 
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Marcosson, ISAAC FREDERICK. Anaconda. 1957. 
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The True Book About Clocks. 1957. 
Modern Clocks, their Repair and Maintenance. Ed. 3, 1955. 


Bruton, Eric. 
ROBINSON, THOMAS RUFER. 


INDUSTRIAL MANAGEMENT 


Production Handbook. Ed. 2. 1958. 
SmiTH, J. SANDFORD. The Management Approach to Electronic Digital Computers. 1957. 


MANUFACTURE 
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PHYSICS 


Louis, PRINCE DE. Physics and Microphysics. 1955. 
BucHHOLz, HERBERT. Elektrische und Magnetische Potentialfelder. 1957. 
Bupd, Acoston. Theoretische Mechanik. 1956. 
GLASSTONE, SAMUEL. Sourcebook on Atomic Energy. Ed. 2. 1958. 
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MeetHaM, A. R. Basic Physics. 1957. 
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YounGER, JoHN Exttiotr. Advanced Dynamics. 1958. 

PLASTICS 


ButzKo, Ropert L. Plastic Sheet Forming. 1958. 
FLoyp, Don EpGar. Polyamide Resins, 1958. 
ROCKETS 


Bates, Davip RoBERT. Space Research and Exploration. 1957. 


Carpin, Martin. Vanguard! 1957. 
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1958. 
SCIENCE 
Conant, JAMES Bryant. Harvard Case Histories in Experimental Science. 1957. 


DO YOU KNOW OF THE FRANKLIN INSTITUTE’S 
NEED FOR GIFTS? 


Since 1824 The Franklin Institute, through wise and conservative management, 
has been partially self-supporting. 

As an old institution, established in the public interest, The Franklin Institute 
has tried to earn most of its cost by diligent and productive work. Like many 
other institutions of a non-profit nature, The Institute has, in these years of an 
expanding economy and rising spiral of costs, found it increasingly difficult to perform 
the broad educational services of which the growing nation and territories are ever 
in need. 

Despite an erroneous public impression gained, perhaps, by the excellence and 
wide range of services performed, that it is a wealthy, richly endowed organization, 
The Institute always needs funds. 

There is a warm satisfaction in giving financial support to an organization that 
will use it constructively to help advance and improve our type of civilization. 

Your gift or bequest, large or small, will be deeply appreciated. We shall be 
pleased to give any additional information regarding gifts or memorials. 

Where property is transferred, title should be in the name of The Franklin 
Institute of the State of Pennsylvania for the Promotion of the Mechanic Arts. 

Please write to: The Secretary, The Franklin Institute, Benjamin Franklin 
Parkway at Twentieth Street, Philadelphia 3, Pennsylvania. 
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THE FRANKLIN INSTITUTE LABORATORIES 
NICOL H. SMITH, DIRECTOR 


EXPERIMENTAL STUDY OF STRESSES 
BY 
ROBERT E. PETERSEN! 


The model used for these tests is scaled from the 30-in. outlet nozzles 
in the APDA Fast Breeder Reactor for the Enrico Fermi Power Plant. 
This reactor is sodium cooled and will be subjected to high temperatures 
but relatively low pressures. The thin walls allowed by low pressure 


requirements make it difficult to draw on industrial experience in piping 
and nozzle design since most previous problems deal with high pressure 


applications. 


Instrumentation and model used in the test. 


1 Senior Staff Engineer. 
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The purpose of this investigation is to determine the relationship 
between the maximum stresses and applied piping reactions. A theo- 
retical analysis of stresses in nozzles and in the adjacent vessel wall is 
extremely difficult in view of the complex geometry. In the case of the 
nozzle considered here the geometry is further complicated by the 
faired knuckle, and the presence of the conical transition of the vessel 
wall just above the nozzle. The toroidal fairing was intended to miti- 
gate thermal shock stresses. 

Longitudinal moments, circumferential moments, and radial thrusts 
with respect to the reactor shell were applied to the nozzle. Stresses 
were determined from strain measurements along the intersections of 
the vessel wall and a horizontal and vertical plane through the axis of 
the pipe. Strain measurements were taken in the longitudinal and cir- 
cumferential directions with respect to the vessel and pipe. These 
strains were measured on both the outer and inner surfaces. From 
these data and from given pipe reactions, we can estimate the stresses 
in the prototype. 

Results of these tests indicate that the nozzle stresses are lower than 
we would have anticipated basing our judgment on Bijlaard’s? analyses. 
There are two unique features of the test model which affected the 
ratios of maximum stress to applied loads; hence the relative roles of 
the vessel transition knuckle and the toroidal fairing of the nozzle in 
decreasing the stresses in the shell cannot be compared. It would be 
desirable to evaluate the toroidal fairing of the nozzle in a vessel with- 
out the conical transition. 


Acknowledgment 


The work described here was done under the sponsorship of Atomic 
Power Development Associates, Inc. 


2 Bijlaard, P.P., “‘Stresses for Radial Loads and External Circumferential and Longitudinal 
Moments in Cylindrical Pressure Vessels’’ ; Report to the Pressure Vessel Research Committee, 
March 1955. 


BOOK REVIEWS 


FEEDBACK CONTROL SysTEMs, by Otto J. M. 
Smith. 694 pages, diagrams, 6 X 9 in. 
New York, McGraw-Hill Book Co., Inc., 
1958. Price, $13.50. 


The sociological phenomenon of the “vi- 
cious circle’ has long been recognized as 
positive feedback. The stabilizing effect of 
negative feedback has been advocated even 
longer as illustrated by the quotation “If 
anyone strikes you on the right cheek, turn 
the other to him, and if anyone sues for your 
shirt, give him your coat too.” To intelli- 
gently apply controls to such feedback sys- 
tems, there is a need for developing general 
methods of closed loop analysis and synthesis. 
This book is written, in the author’s words, 
as ‘“‘a unified philosophy for the analysis and 
design of all types of feedback systems.” He 
includes social, economic, government, and 
biological systems with the more clearly de- 
fined physical ones. 


The book is designed as a reference and 
textbook for graduate study, assuming prior 
courses in basic feedback control and statis- 


tical analysis. The four main divisions of the 
text are: linear analysis, linear synthesis, 
non-linear analysis and non-linear synthesis 
and each is treated extensively. The first 
quarter on linear analysis begins with the 
conventional treatment of a single feedback 
loop. It expands into the analysis of multiple 
loops and a particularly lucid description of 
power spectra. 

The second section, dealing with linear 
synthesis, opens with the definition of the 
author’s choice of an optimum system, 
namely, one with minimum error power. The 
importance of a statistical description of the 
expected input is stressed. A new function 
is introduced, called the “realizability opera- 
tor’’ which removes, somewhat magically, all 
unrealizable poles in the optimum design. 
Coming back down to earth for an interval, 
the author presents an excellent section on 
design techniques, indicating an intimate 
knowledge of the mechanics of synthesis by 
graphical methods. These techniques are 
extended to cover the simpler non-linear sys- 
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tems having dead time and distributed 
parameters. 

The third and fourth sections present what 
is probably the most comprehensive approach 
to non-linear systems to date, even consider- 
ing the excellent work of John Truxal. Relay 
systems with and without multilevel signals 
are treated at length by the method of two 
loci in the L-plane. Chapters on unilateral 
and bilateral non-linearities and discontinui- 
ties follow. Examples of these are saturation 
hysteresis and dead zone as the unilateral 
type, and backlash and coulomb friction as 
the bilateral type. Both may lead to dis- 
continuities such as jump phenomena and 
lock-up oscillation. 

Non-linear synthesis is handled primarily 
by the phase-plane method, graphically and 
with an eye toward real components. Pre- 
dictor control for compensating non-linearities 
is described and a general treatment of carrier 
systems completes the section. The appendix 
contains a set of compensating networks with 
their mathematical and graphical description. 

The book is an outstanding contribution to 
the field, being valuable as a text for advanced 
graduate study as well as a working reference 
for the designer of complex systems. It 
effectively fills out the McGraw-Hill Series in 
Control Systems Engineering. 

CHARLES A, BELSTERLING 
The Franklin Institute Laboratories 


DISLOCATIONS AND MECHANICAL PROPERTIES 
or CrystALs, edited by J. C. Fisher, W. G. 
Johnson, R. Thomson and T. Vreeland. 
634 pages, illustrations, 6 X9 in. New 
York, John Wiley & Sons, Inc., 1957, 
Price, $15.00. 

This book is the result of a small inter- 
national conference held Sept. 6-8, 1956 at 
Lake Placid, New York initiated by N. F. 
Mott and J. H. Hollomon who recognized the 
need to review the present state of research 
both theoretical and experimental on disloca- 
tions and their relationship to the mechanical 
behavior of crystals. It represents the final 
results of considerable recent work on dislo- 
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cations in crystals contributed by fifty-four 
scientists active in this field. Individual 
papers and summaries of papers which were 
presented at the conference are included in 
the book together with edited versions of the 
discussions which followed each paper. The 
discussions contributed much in pointing out 
the need for further work in this field due toa 
lack of agreement among the participants and 
complete understanding of many of the associ- 
ated phenomena, 

The book has been logically subdivided into 
eight chapters which include the following: 
1, Direct Observations of Dislocations; II, 
Deformation of Pure Single Crystals; III, 
Work Hardening and Recovery; IV, Alloy 
Crystals, Impurities, Yield Point Phenomena; 
V, Dislocation Damping and Fatigue; VI, 
Theory of Dislocations; VII, Whiskers and 
Thin Crystals; and VIII, Radiation Damage. 

Considerable work in the observation of 
dislocations in transparent KCI crystals was 
contributed by S. Amelinckx. Numerous 
photographs showing dislocations and the 
pinning of dislocations into nets are presented 
and analysed. The results show in a very 
striking manner the pinning of dislocations by 
precipitates in addition to the intersection, 
annihilation and climb of dislocations. Other 
papers describing techniques for observing 
dislocations in Silicon, Silver Halides, Lithium 
Fluoride, Aluminum, Silver and Iron-nickel 
alloys are given. Direct observations of dis- 
location markings were generally supported 
by studies of their movement under applied 
stresses. 

The subject of deformation of face-centered 
cubic metals by twinning has been a contro- 
versial one in recent years and a very inter- 
esting paper by T. H. Blewitt, R. R. Coltman 
and J. K. Redman on the deformation of 
copper single crystals at low temperatures was 
presented in which twinning was observed in 
crystals pulled in tension at 77.3° K. Other 
work of the deformation of pure single crystals 
involved the subject of surface sources and 
their importance in the onset of plastic flow. 

A very detailed and comprehensive analysis 
of the mechanism of glide and work hardening 
in face-centered cubic and hexagonal close 
packed metals based upon a dislocation model 
is presented by Alfred Seeger in which the 
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three stages easy glide, work hardening and 
recovery are discussed. 

The theory of hardening due to the elastic 
and chemical interaction of solute atoms with 
dislocations and stacking faults is often used 
to explain the yield strength of binary alloys. 
The work of Hideji Suzuki presented in Chap- 
ter IV deals with the evaluation of the yield 
strength under various conditions of elastic 
and chemical interaction and compares them 
with experimental observations. 

Internal friction phenomena and fatigue 
have also been analyzed in terms of a disloca- 
tion model but due to the complexity of these 
processes advances in this field have been im- 
peded. Theoretical models capable of de- 
scribing some aspects of damping phenomena 
have been proposed. The behavior of metals 
under reversed stresses is treated by N. F. 
Mott in which an attempt is made to improve 
current theory of work hardening under cyclic 
stress. 

Additional important contributions to the 
theory of dislocations are made by G. Lieb- 
fried on the thermal motion of dislocations, 
by Alfred Seeger on the nonlinearity of the 
elastic laws in the vicinity of a dislocation, by 
A. H. Cottrell on the intersection of gliding 
screw dislocations and by J. C. Fisher on 
prismatic loops as Frank-Reed sources. Dis- 
cussions of these papers by other noted 
authorities in the field proves to be extremely 
interesting and enlightening. 

Radiation damage in metals is the subject 
of the final chapter in which the study of in- 
terstitials and vacancies is discussed. Al- 
though far from complete some of the more 
important aspects of cavity formation and its 
effect on electrical and mechanical behavior 
is covered. 

Many other papers are presented and dis- 
cussed in the book which will be of consider- 
able value to many workers active in the field. 
This voiume is not a text book from which to 
study the complete theory of dislocations. It 
does, however, represent the accumulation of 
experimental work, analyses and thinking of 
numerous scientists presently active in the 
field of dislocations and will best serve the 
researcher or scholar already well versed and 
familiar with the subject of dislocations. 


Victor DAMIANO 
The Franklin Institute Laboratories 
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THREE DIMENSIONAL Dynamics, by C. E. 
Easthope. 277 pages, diagrams, 5$ X 84 
in. New York, Academic Press, Inc., 1958. 
Price, $7.80. 


Although there are still colleges in this 
country where a student can graduate in 
physics without ever coming into serious con- 
tact with vector analysis, they are definitely 
the exception instead of the rule. Even in 
engineering schools, vector treatment of elec- 
tromagnetics is common, and that of dynamics 
becoming so. 

A few years ago there were almost no suit- 
able texts which joined vector analysis to tech- 
nical mechanics. There still are not many, 
but they are appearing at an increasing rate. 
Unfortunately, the present work, despite its 
promising title, will probably find little use in 
American schools. An intermediate text, it 
covers material up to Euler’s equations but 
not beyond. The exposition is entirely 
mathematical, after the British manner. Al- 
though the author’s language is lucid and 
relaxed, his excessive formality, highly arti- 
ficial examples, meager illustrations, and 
proof-type problems should prove grave 
handicaps in a U. S. classroom. Perhaps 
this text would be most useful in teaching 
vector methods to students already familiar 
with mechanics. 

C, YAGER 
The Franklin Institute Laboratories 


ELECTRON TUBES AND SEMICONDUCTORS, by 
Joseph J. DeFrance. 288 pages, illustra- 
tions, 6 X 9 in. Englewood Cliffs, N. J., 
Prentice-Hall, Inc., 1958. Price, $9.00. 


This is a very elementary book serving to 
introduce tubes and semiconductor devices to 
persons who have little acquaintance with the 
engineering aspects of radio or electronics. A 
similar book was published in 1936 by E. D. 
McArthur, but this earlier work contained 
more fundamental information relating to 
physical principles. In that book Richard- 
son’s equation was presented and the improve- 
ment made upon it by Dushman, so that the 
student obtained some solid background of 
electronic principles. He was thereafter able 
to discuss these effects by name and show 
some intelligence when, for example, the 
Schottky effect was mentioned. 
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Joseph DeFrance’s book on the other hand 
appears to be primarily a collection of tube 
manual illustrations and characteristic curves. 
The explanations also bear a resemblance to 
those presented in the opening pages of the 
tube manuals or advertising literature. Most 
technicians are surprisingly familiar with this 
material and usually possess a large degree of 
knowledge concerning the practical aspects of 
the associated circuits. 

Since the book is aimed at the ‘technical 
institute’ or technician level, it would seem 
better to teach some basic electronic prin- 
ciples, for it is here that the technician’s know- 
ledge is most frequently lacking. On the 
other hand, the kind of information which this 
author has chosen for his book is usually 
acquired by technicians from informal sources. 
The volume might bridge certain gaps in the 
technician’s general knowledge, for it does 
cover in a comprehensive manner the newer 
subjects such as transistors and recent ad- 
vances in photoelectric devices. However, it 
is difficult to see how such qualitative infor- 
mation could substantially increase one’s 
effectiveness at the working level. 


C. W. HARGENS 
The Franklin Institute Laboratories 


ELEcTRICAL DISCHARGES IN GASES, by F. M. 
Penning. 75 pages, diagrams, 6 X 9 in. 
New York, The Macmillan Co., 1958. 
Price, $3.00. 


In the introduction, W. de Groot explains 
that this is a translation of a book published 
in Dutch in 1955, two years after the author’s 
death. He says further that it offers “a con- 
cise synopsis which can profitably beused . . . 
as a basis for further study. The contents 
naturally reflect to a great extent the nature 
and trend of the investigations carried out in 
the Philips Eindhoven Laboratories, in which 
investigations Dr. Penning had such a very 
important share.” 

This little book (including index and bibli- 
ography, it is only 75 pages long) explains, in 
descriptive terms, the underlying physical 
mechanisms of different kinds of direct current 
discharges simply and clearly. Not a com- 
prehensive treatise, it makes no pretense of 
covering every possible topic, nor does it pre- 
sent any complicated mathematics. Never- 
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theless, it contains a wealth of detail, and 
includes equations and numerical data to a 
sufficient extent to give the reader a feel for 
the magnitude of the important quantities. 
Even the reader who has only a very sketchy 
notion of the nature of electrical discharges 
will find the explanations easy to follow, con- 
ceptually illuminating, and presented with 
stylistic and pedagogical excellence. 


G. P. WACHTELL 
The Franklin Institute Laboratories 


ADVANCES IN ELECTRONICS AND ELECTRONIC 
Puysics, Vol. IX, edited by L. Marton. 
347 pages, illustrations, 6 X9 in. New 
York, Academic Press, Inc., 1957. Price, 
$9.00. 


The advent of the International Geophysi- 
cal Year has sparked research in many fields 
but in perhaps no one are the advances so 
significant as in the field of electronics. Not 
only have advances been achieved in instru- 
mentation application and the use of electron- 
ics as a tool, but also as probes in analyzing 
terrestrial and atmospheric phenomena. 

Five contributors cover subjects dealing 
with the atmosphere, the Aurora Borealis, 
Negative Ions, Radio Observations of Meteors, 
Intensity Variations in Cosmic Rays and a 
review of Radio-Wave Propogation comprise 
this section. The remaining two sections deal 
with Electronics in Oceanography and Con- 
tributions of Electronics to Seismology and 
Geomagnetism. 

The variety of topics in a single volume 
makes it difficult to assess this book. This 
reviewer has had professional interest in the 
aurora, meteors and seismology. Thus, out 
of seven sections, only three are of sufficient 
importance to maintain interest. If the book 
were intended primarily for electronic en- 
gineers and physicists there isn’t enough meat 
in these sections. It is unfortunate that Dr. 
Marton tried this experiment of using a theme 
and crossing disciplines instead of restricting 
the book to a single discipline. 

Despite these shortcomings Advances in 
Electronics and Electronic Physics contains 
some excellent presentations and is a most 
welcome addition to our libraries. 


I. M. Levitt 
Fels Planetarium of The Franklin Institute 
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OrGANIC ELECTRODE PRocEssEs, by Milton 
J. Allen. 174 pages, diagrams, 53 X 8} in. 
New York, Reinhold Publishing Corp. 
1958. Price, $6.50. 


The reviewer cannot help but compare this 
book to the section on “Electrolytic Reac- 
tions’ by Sherlocke Swann, Jr. in “Technique 
of Organic Chemistry” Volume II, second 
edition (Interscience, 1956). Dr. Swann’s 
treatment contains an extensive table of elec- 
trochemical reactions which makes the data 
easily accessible to the reader with a minimum 
of text. He includes a number of illustrative 
procedures for typical reactions and a brief 
discussion of electrochemical theory. Dr. 
Allen’s book presents essentially the same in- 
formation less economically, but does contain 
more practical information on instrumenta- 
tion and technique of electrolysis. 


M. M. LaBes 
The Franklin Institute Laboratories 


VIBRATION AND Impact, by Ralph Burton. 
310 pages, diagrams, 6 X 9 in. Reading 
(Mass.), Addison-Wesley Publishing Co., 
Inc., 1958. Price, $8.50. 


In this introductory text, written for 
advanced undergraduates, the author has 
intended to convey an understanding of oscil- 
latory phenomena in general, as an introduc- 
tion to the more specialized subjects of 
interest to mechanical engineering students. 
He succeeds in making his presentation attrac- 
tive through a judicious combination of theory 
and practical applications. 

The history of vibration study, and a num- 
ber of examples illustrating its importance, 
are briefly outlined in the first chapter. Suc- 
ceeding chapters deal with the basic topics of 
free vibrations of systems having one degree 
of freedom; typical vibration problems en- 
countered in machinery ; damping; and steady 
forced vibrations. The analysis of periodic 
non-harmonic forcing functions leads over into 
a discussion of impact and other forms of 
transient loading, by way of Fourier series. 
Nonlinear vibrations are discussed in Chapter 
7, and systems with two degrees of freedom 
are dealt with in Chapter D, which contains 
a brief but adequate demonstration of the 
orthogonality of the principal modes. 

Chapter 8 presents a useful digression into 
vibration measuring instruments and the 
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simpler electrical analogs of mechanical sys- 
tems. A related chapter, No. B, deals briefly 
with control and feedback systems. The 
analytical methods employed here are admit- 
tedly limited, and do not allow the ease of 
manipulation permitted by the more standard 
Laplace transform methods. 

Two standard methods of computing 
natural frequencies of systems with many 
degrees of freedom are described in Chapter 
10. Analyses of vibrations in continuous 
systems (including gas columns, strings, 
beams and plates) are outlined. A final 
chapter discusses the relationships between 
vibration and fatigue phenomena. 

A summary, list of references and problems 
(with solutions) conclude each chapter. 
Many problems are worked out in full, in 
keeping with the pedagogic intent of the book. 
Much trouble has been taken with the illus- 
trations, although confusion sometimes results 
from a combination of isometric and diagram- 
matic views. On the whole, however, this 
book, on account of its generally clear pres- 
entation and interpretation of physical 
phenomena, can be highly recommended as 
an introduction to an increasingly important 
subject. 

Kiaus CAPPEL 
The Franklin Institute Laboratories 


THE PENTAERYTHRITOLS, by Evelyn Berlow, 
Robert H. Barth and John E. Snow. ACS 
Monograph 136, 317 pages,6 X 9in. New 
York, Reinhold Publishing Corp., 1958. 
Price, $10.00. 

The brief history of pentaerythritol might 
read, ‘‘white crystal—discovered 1882; white 
powder—multi million dollar industry, 1958.” 
The properties, compounds, and uses of this 
material are summarized and referenced in 
this volume. 

The wide number of uses for pentaerythritol 
stems from the high reactivity of its four 
hydroxyl groups. Thus, many compounds 
are available for such diverse tasks as filling 
charges for artillery shells, plasticizers to im- 
prove the properties of plastics, preparation 
of synthetic drying oils for paints and var- 
nishes, and a base material for insecticides. 
Further uses for fatty acid esters of the 
material may be in high temperature lubri- 
cants for turbojet and turboprop engines. 
Discussions of preparations and properties 
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of the compounds of the pentaerythritols are 
presented according to type. Each chapter 
contains a catalog of compounds and is well 
referenced with as many as nine pages of 
bibliography in the case of nitrates. It is 
somewhat upsetting, however, to find the uses 
for pentaerythritol scattered through the 
chapters. Such format will, of necessity, 
force the reader who is interested in applica- 
tions of the material to ferret out the informa- 
tion. It would seem that a separate chapter, 
inserted possibly after the chapter on prepara- 
tion of the parent compound, would be 
desirable. 

The Pentaerythritols is not a book for the 
reader of general chemical interest. It should, 
however, be recommended for the chemist 
engaged in work involving the use or manu- 
facture of the material or its compounds, and 
it will undoubtedly be a valuable addition to 
the shelves of the technical reference library. 


RICHARD HOLLINGER 
The Franklin Institute Laboratories 


HEAT TRANSFER, Volume II, by Max Jakob 
with the technical and editorial assistance 
of Stothe P. Kezios. 652 pages, illustra- 
tions, 6 X 9 in. New York, John Wiley & 
Sons, Inc., 1957. Price, $15.00. 


The multitude of devoted students, col- 
leagues, and friends of Dr. Max Jakob who 
mourned his death in January 1955 were 
saddened additionally by the prospect that 
the much awaited second volume of Heat 
Transfer would never be published. The 
scientific community, therefore, owes a debt 
of gratitude to Dr. Jakob’s wife Anna who had 
the desire to have the volume made ready for 
publication and especially to Dr. Kezios whose 
efforts made that desire a reality. Both vol- 
umes stand asa lasting contribution by a kind, 
humble, and learned man. 

Volume II follows the same formal organi- 
zation presented in Volume I with the sec- 
tions, chapters, and equations being numbered 
consecutively through both volumes. The 
table of contents contains the many subdivi- 
sions comprising each chapter so that a specific 
development within a major topic may be 
located quickly and easily. Of significant 
value to the student is the detail paid to the 
historical aspect of the development of partic- 
ular topics. 
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The first two chapters of the book deal with 
thermal radiation through absorbing and non- 
absorbing media. In Volume I the laws of 
thermal radiation and the influence of material 
properties on radiation were presented with- 
out considering the medium between two 
radiating surfaces; geometrical factors were 
mentioned but not analyzed. In Volume II 
the medium on and between two radiating 
surfaces as well as the geometry of the surfaces 
is treated in a complete manner. Virtually 
all possible conditions involving simple con- 
figurations are analyzed in a manner suffi- 
ciently rigorous to satisfy mathematicians, 
physicists, and engineers. Certain recent 
developments, for example, the determination 
of heat transfer from strongly shocked air are 
not discussed but undoubtedly will be in- 
cluded in later editions or supplements. Ra- 
diation to and from gases is restricted to sys- 
tems comprising carbon dioxide and/or water 
vapor. 

The major portion of the book has been 
devoted to selected fields of application. 
Heat transfer in thermometry, heat ex- 
changers, regenerators, boundary layer heat 
transfer, and heat transfer in liquid metals are 
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treated in some detail. Of special interest to 
those working with high fluid velocities are 
Chapters 38, 39 and 40 which deal with heat 
transfer through laminar and _ turbulent 
boundary layers as well as the fundamentals of 
transpiration cooling. A thorough discussion 
of analytical and experimental work up to 
1953 is presented. 

The last part of the book comprises supple- 
ments to various chapters of Volume I in- 
cluding some recent data on physical proper- 
ties of liquid metals. Additional problems 
relating to subject materials in Volume I are 
also included. 

An extensive bibliography and author index 
is presented for the material covered in the 
book. The researcher’s task is made much 
simpler by virtue of the completeness of the 
reference list. 

This reviewer recommends the book, indeed 
both volumes, for the engineer or scientist who 
requires a reference to the fundamental laws 
and theories of heat transfer along with suffi- 
cient applications to indicate the validity and 
limitations of the theories. 


Exuiott L. Katz 
The Martin Company 


ANALYTICAL MECHANICS, by F. W. Van Name, 
Jr. 171 pages, diagrams, 6 X 9 in. New 
York, Prentice-Hall, Inc., 1958. Price, 
$5.75. 


A readable, elementary account of dy- 
namics for physics majors, presupposing no 
more than general physics and calculus. De- 
signed expressly as a teaching text, it includes 
only material which can be effectively covered 
inonesemester. Vectorsare introduced at the 
beginning and used throughout. 


Atomic PHysics AND HUMAN KNOWLEDGE, 
by Niels Bohr. 101 pages, 5} X 9 in. 
New York, John Wiley & Sons, Inc., 1958. 
Price, $3.95. 


This is a collection of seven articles by the 
Nobel prize winner whose earlier work forms 
the basis of present atomic physics. Four of 
the articles have been published during the 
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period from 1933 to 1955; the other three are 
addresses delivered on various occasions but 
not published previously. 

These articles deal with the impact which 
atomic physics has had on human knowledge, 
in such fields as biology, anthropology, psy- 
chology, natural philosphy and physiology. 
The volume, representing the attitude of one 
of this century’s greatest scientific minds, will 
be welcomed as a lucid treatment of a subject 
which profoundly affects us all. 


Wauat’s HAPPENED TO Our HiGH ScHOooLs?, 
by John Francis Latimer. 196 pages, 6 x 9 
in. Washington, D. C., Public Affairs 
Press, 1958. Price, $3.25. 


Appropriately published at a time when the 
United States is taking shocked inventory of 
the deficiencies in its educational system, this 
book has been several years in preparation. 
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The author, Assistant Dean of Faculties at 
George Washington University, has gathered 
the statistics for his book from many sources, 
including the U. S. Office of Education and 
the National Education Association. 

The book opens with a history of public 
education in the United States and proceeds 
to give data on enrollment and subjects taught 
for the years 1900, 1910, 1922, 1934 and 1949, 
showing the extension of educational oppor- 
tunities and expansion of the curriculum. 

In a chapter called “Split-Level Educa- 
tion,”’ the author discusses the division of 
secondary education into two groups—one for 
college preparatory students and one for those 
not planning to attend college. The gradual 
development and separation of the two cur- 
ricula are described in detail, with sections 
devoted to the Dewey philosophy that the 
aim of education should not be exclusively 
intellectual, but social as well. 

The effect of this split-level curriculum is 
given in Chapter VI, in which the author 
traces the decline in enrollment for such sub- 
jects as mathematics and the sciences. This 
chapter also carries the study into the college 
level, with tables of majors listed for the 
period 1948-1955, reflecting the trend away 
from the “‘classical’”’ subjects. 

In the final chapter, the author appraises 
the situation and offers several suggestions, 
among them one that the students be grouped 
according to ability, but that the less gifted 
child’s education should differ only in quan- 
tity, not in kind, from that of the gifted child. 

This is a solidly documented work, on a 
basic problem facing us today. It should be 
of intense interest to teachers, educationists, 
parents, school board members and faculty 
members of colleges and universities. 


APPLIED MATHEMATICS FOR ENGINEERS AND 
Puysicists, by Louis A. Pipes. Second 
edition, 723 pages, diagrams, 6 X 9 in. 
New York, McGraw-Hill Book Co., Inc., 
1958. Price, $8.75. 


Substantial additions have been made in 
this new edition of a well known work. It 
remains primarily a reference, but is suitable 
also as a text in the standard omnibus course 
in applied mathematics. Coverage is strong 
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in most topics, notably matrices, but sketchy 
on a few. A substantial bibliography is in- 
cluded, partially annotated. There are many 
problems, but no answers. 


MoperN ComputiING METHops, by members 
of the staff of the Imperial College of Sci- 
ence and Technology. 129 pages, dia- 
grams, 6 X 94 in. New York, Philosophi- 
cal Library, 1958. Price, $8.75. 


This text renders a condensed account of 
general computing methods with emphasis 
being placed on practical and simple problems 
and worked examples. Algebraic equations 
and matrices, finite difference methods and 
ordinary and partial differential equations are 
topics dealt with briefly in the twelve chapters. 
Included also is an extensive bibliography. 


THERMODYNAMICS OF HEAT Power, by Virgil 
Moring Faires. Second edition, 432 
pages, illustrations, 6 X 9}in. New York, 
The Macmillan Co., 1958. Price, $8.00. 


This is a second edition, appreciably re- 
vised, of the author’s Theory and Practice of 
Heat Engines. Intended as primarily a 
theory book, it nevertheless contains much 
technological detail, including material on 
heat exchange in reactors. Following the 
pattern of Professor Faires’ other very teach- 
able texts, there are plenty of excellent illus- 
trations and a wealth of numerical problems, 
almost all with answers. 


TECHNIQUES OF PRESSWORKING SHEET 
METAL, by Donald F. Eary and Edward A. 
Reed. 472 pages, diagrams, illustrations, 
6X9 in. New York, Prentice-Hall, Inc., 
1958. Price, $12.00. 


A reference book dealing thoroughly with 
the processes of sheet metal work. Based on 
actual plant experience, this text outlines and 
defines for clarity the design and construction 
of dies. The authors analyze vector forces 
and use strength of materials to determine die 
component size, with a high degree of detail 
and an extensive coverage of terminology. A 
vast number of photographs, line drawings, 
charts and tables illustrate clearly much of 
the written text. 
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PUBLICATIONS RECEIVED 


The following books have been received in the JouRNAL office for review. Lack of space 
prevents publication of more than this listing. 


FLurp MEcuanics, by Victor L. Streeter. Second edition, 480 pages, diagrams, 6 X 9 in. 
New York, McGraw-Hill Book Co., 1958. Price, $7.50. 

CHEMICAL PUBLICATIONS, THEIR NATURE AND Usk, by M. G. Mellon. Third edition, 327 
pages, illustrations, 6 X 9 in. New York, McGraw-Hill Book Co., 1958. Price, $7.00. 

CoLLEGE CHEMISTRY, by Paul R. Frey. Second edition, 707 pages, illustrations, diagrams, 
6 X9in. New Jersey, Prentice-Hall Inc., 1958. Price, $6.95. 

Tue GOLDEN NuMBER, by M. Borissavlievitch. 91 pages, diagrams, 54 X 7 in. New York, 
Philosophical Library, 1958. Price, $4.75. 

A SHort History or ANATOMY & PHYSIOLOGY FROM THE GREEKS TO HARVEY, by Charles 
Singer. Second edition, 209 pages, diagrams, 54 X 8 in. New York, Dover Publications 
Inc., 1957. Price, $1.75 (paper). 

Tue Story OF ARCHAEOLOGY, by Agnes Allen. 245 pages, diagrams, 5} X 7}in. New York, 
Philosophical Library, 1958. Price, $4.75. 

ARTIFICIAL STIMULATION OF RAIN, edited by Helmut Weickmann and Waldo Smith. 427 
pages, illustrations, diagrams, 74 X 9} in. New York, Pergamon Press, 1957. Price, 
$15.00. 

MECHANICAL ENGINEERING LABoratory, by Charles W. Messersmith, Cecil F. Warner and 
Robert A. Olsen. Second edition, 179 pages, illustrations, diagrams, 8 X 11 in. New 
York, John Wiley & Sons, Inc., 1958. Price, $5.50 (paper). 

THE ILLUMINATED Book, Its History AND Propuction, by David Diringer. 524 pages, 
plates, illustrations, 6 x 93 in. New York, Philosophical Library, n.d. Price, $25.00. 

THE DIRECTION OF RESEARCH ESTABLISHMENTS, a symposium. v.p., illustrations, 64 X 9 in. 
New York, Philosophical Library, 1957. Price, $12.00. 

Puysics PROBLEMS AND How To SoLvE THEM, by Clarence E. Bennett. 240 pages, diagrams, 
54 X 8} in. New York, Barnes & Noble, Inc., 1958. Price, $1.75 (paper). 

ELEMENTARY PRACTICAL ORGANIC CHEMISTRY, Part 3 Quantitative Organic Analysis, by Arthur 
I. Vogel. Pages 645 to 840, diagrams, 54 X 8} in. New York, Longmans, Green & Co., 
Inc., 1958. Price, $4.50. 

CLIMATOLOGICAL SUMMARY MouLp Bay, N. W. T., CANADA. 92 pages, diagrams, 8 X 11 in. 
Canada, Department of Transport, Meteorological Division, 1957. No price, (paper). 

PHILOSOPHY OF MATHEMATICS AND NATURAL SCIENCE, by Hermann Weyl. Third printing. 
311 pages, diagrams, 6 X 9 in. Princeton, Princeton University Press, 1958. Price, 
$1.95 (paper). 


hy 
in- 
rs 
a- 
of 
is 
1S 
is 
d 
il 
2 
f 
a 
1 
1 
: 


CURRENT TOPICS 


Rockets Probe Upper Atmosphere 
for IGY.—A highly successful season 
of rocket probings of the upper atmos- 
phere was completed in April at Fort 
Churchill, Manitoba, Canada, as part 
of the US International Geophysical 
Year program. Reports reaching the 
National Academy of Sciences’ IGY 
Committee in Washington indicate 
that 22 of the 23 firings in 1958 pro- 
duced the desired scientific data. 
Since IGY began in July 1957, 41 
rockets have been launched at Fort 
Churchill, the most extensive rocket 
sounding program ever conducted in 
the higher latitudes. Firings there 
were resumed in mid-June, with at 
least 32 more scheduled before IGY 
ends. 

The rocket firings explored the area 
between 46 and 160 miles above the 
earth’s surface, a portion of the atmos- 
phere which is out of reach for bal- 
loons and, for the most part, below 
the level of satellite orbits. Eleven 
types of experiments were conducted 
in this ionosphere layer, which, in ad- 
dition to being the subject of basic 
scientific research, is of immediate 
practical interest because of its effect 
on long distance radio communica- 
tions. Included were the first suc- 
cessful firings into visible auroras. 

Data now being processed are ex- 
pected to throw new light on such 
matters as the altitude of the electric 
currents responsible for disturbances 
in the earth’s magnetic field during an 
aurora and on the nature, extent, and 
location of the electrons and ions at 
high altitudes, which are governing 
factors in radio transmission. 

Ingenious techniques were used to 
measure the vertical distribution of 
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pressure, temperature and density of 
this segment of the atmosphere. 
Paradoxially, these techniques were 
made necessary by the rocket itself. 
Although the rocket is the only means 
of attaining the right heights, its great 
speed does not allow an ordinary ther- 
mometer probe to come into equilib- 
rium with the atmosphere. It there- 
fore becomes necessary to calculate 
temperature indirectly from pressure, 
density, or velocity of sound. 

In the sound-velocity method, 18 
explosive grenades are ejected from 
the rocket at timed intervals. On the 
ground, an array of sensitive micro- 
phones is spaced 1000 yards apart in 
the form of across. Each microphone 
“hears’’ the burst at a slightly differ- 
ent time. By noting the time taken 
for the sound to travel from rocket to 
ground, average temperature of the 
air can be computed, and the differ- 
ence in arrival time at the various 
microphones yields a measure of av- 
erage upper level wind velocity and 
direction. 

The Fort Churchill operation is the 
result of the combined efforts of many 
agencies. The facility itself was es- 
tablished jointly by the Canadian and 
US National Committees for the IGY 
and the defense establishments of both 
countries. The U.S. Army is respon- 
sible for the complicated ground sup- 
port phase—firing of just two rockets 
involved 4,408 man-hours of overtime 
in activities ranging from road blocks 
for safety to actual launchings, for 
which the Navy provides firing crews. 

Scientific participants from the 
U. S. include the Air Force Cam- 
bridge Research Center, Ballistics 
Research Laboratory of Aberdeen, 
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Md., Signal Engineering Laboratory 
of Fort Monmouth, N. J., State Uni- 
versity of lowa, U.S. Naval Research 
Laboratory, which with assistance of 
personnel from New Mexico College 
of Agriculture and Mechanic Arts op- 
erates an ionosphere station and tele- 
metering trailers, the University of 
Michigan, and the University of Utah. 
Canadian organizations make auroral 
observations, magnetic recordings and 
ionosphere soundings. 

Rocket vehicles used or planned are 
the Loki Phase II-Dart, the Nike- 
Cajun and Nike-Asp, and the Aerobee 
and Aerobee-Hi. The Aerobees are 
the largest, and carry 150 pounds of 
scientific instrumentation. 

Other phases of the US-IGY rocket 
program are carried on at White Sands 
and Alamogordo, N. Mex., Pt. Mugu, 
Calif., and from shipboard in the 
Arctic, equatorial regions and the 
Antarctic, where the rockets are 
launched from balloons. 


Electronic Highway of Tomorrow in 
Operation.—A segment of the “elec- 
tronic highway of tomorrow’ went 
into operation in May at the David 
Sarnoff Research Center of the Radio 
Corporation of America at Princeton, 
N. J. Cars entering the grounds of 
the research center, headquarters of 
RCA Laboratories, passed over elec- 
tronic loops buried in the entrance 
road. Transistorized detector units 
connected to the loops counted each 
of the passing cars, measured the 
speed of each, and flashed a polite 
warning saying ‘‘Slower Please’’ to 
every driver exceeding the posted 
speed limit. 

The new equipment is a transistor- 
ized variation of the comprehensive 
electronic vehicle control system de- 
veloped earlier by RCA and demon- 
strated last year at Lincoln, Nebraska, 
in cooperation with the Nebraska De- 
partment of Roads. At Lincoln, the 


CuRRENT Topics 75 


system showed its ability to guide 
cars automatically along the highway 
and to provide electrical signals for 
the automatic control of braking and 
steering. 

The new RCA Laboratories instal- 
lation is designed to demonstrate 
further immediate applications of the 
system in specific traffic control jobs. 
Dr. Vladimir K. Zworykin, Honorary 
Vice-President of RCA, directed de- 
velopment of the entire vehicle con- 
trol system, in cooperation with an 
RCA Laboratories research team in- 
cluding L. E. Flory, W. S. Pike, and 
G. W. Gray. 

The distinction of receiving the 
first speed warning fell to Mr. Flory, 
who purposely drove past the 20- 
miles-per-hour speed limit sign at the 
research center entrance at a smart 
30-miles-per-hour clip. As the car 
passed over the two wire loops buried 
in the pavement, four things happened 
simultaneously : 


1. The transistorized detector units 
connected with the buried loops de- 
tected the passage of the car and 
measured its speed. 

2. Since the speed was greater than 
the specified limit, the circuits auto- 
matically illuminated the warning 
sign at the roadside some twenty 
yards ahead of the car. 

3. A simple counter connected to 
the detector units recorded the pas- 
sage of the car. 

4. A second counter automatically 
recorded the passage of a car in excess 
of the speed limit. 


‘The idea of such an installation on 
the highway would be to study driver 
habits and to determine whether or 
not a posted speed limit is realistic,”’ 
said Dr. Zworykin. ‘By installing 
such a system at the entrance to the 
David Sarnoff Research Center, we 
will have an opportunity to study the 
functioning of the equipment over ex- 
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tended periods and under all condi- 
tions of weather.” 

Ultimately, according to Dr. Zwor- 
ykin, another loop can be installed 
farther along the road to determine 
how many drivers slow down after 
receiving the warning from the auto- 
matic sign, ‘“‘thereby giving us another 
indication of driver reactions.” 

Dr. Zworykin emphasized that the 
newly installed units are components 
that would be standard in the com- 
prehensive electronic control system 
designed for the ultimate automatic 
control of vehicles on the highway. 
Like the demonstration equipment at 
the laboratory entrance, the larger 
system consists basically of wire loops 
buried in the pavement along a high- 
way, with detector units associated 
with each loop to pick up the signals 
generated by passage of cars over the 
loop. In the complete system, these 
signals are fed back through antennas 
under the road to activate automatic 
controls in following cars, or they are 
used to operate various roadside warn- 
ing devices. 

“In the complete system,”’ said Dr. 
Zworykin, ‘“‘any segment at any point 
on the highway can be used to per- 
form the counting and speed meas- 
urement functions demonstrated in 
this laboratory installation. 

“Various radar and mechanical 
techniques are now used for this pur- 
pose, but these existing devices re- 
quire special installation wherever 
they are to be used, and they bear no 
relation to the complete system in the 
highway itself. By using elements of 
the complete RCA system for these 
specific and limited functions, a va- 
riety of installations can be made in 
the near future with the assurance 
that they will function also as integral 
parts of the more extensive highway 
installations of tomorrow.” 

Mr. Flory, as a member of the re- 
search team responsible for develop- 
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ment of the electronic highway con- 
trol system, gave this brief descrip- 
tion of the way in which the newly 
installed RCA Laboratories equip- 
ment operates: 

The detector loops in the roadway 
itself consist of rectangles of wire 
6 X 20 ft. in outline. In operation, 
a voltage from a high-frequency power 
line is applied to each of the loops. A 
vehicle passing over the loop causes 
a variation in the current passing 
through the wire. This variation is 
immediately detected by a _ small 
transistorized detector circuit placed 
at the side of the road. The detector 
circuit then produces an output signal 
which is used to perform the desired 
action. In the RCA Laboratories in- 
stallation, the signals from the detec- 
tor units activate the counters. In 
addition, the time interval between 
the variations in the two loops is 
measured to determine the speed of 
the vehicle, and the illuminated sign 
is turned on if the speed exceeds a pre- 
determined limit. 


Consumable Weld Insert for Thick 
Pressure Vessel Wall.—A consuma- 
ble weld insert, developed originally 
for butt welding pipe, has been used 
successfully by Newport News Ship- 
building & Drydock Co., Newport 
News, Va., to join the hemispherical 
halves of a nuclear reactor pressure 
vessel. 

The pressure vessel presented a 
difficult welding problem in that all 
welding had to be done from the out- 
side while at the same time the inside 
of the weld had to be smooth, uniform 
and crevice free. The 60-in. I.D. 
pressure vessel is a sphere of 4-in. 
mild steel internally clad with an 
0.4-in. layer of type 347 stainless steel. 
In operation, the vessel handles a 
highly corrosive process fluid of heavy 
water under 2000 psi. pressure. 
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Although the use of the weld insert 
requires more careful preparation of 
the joint than standard welding prac- 
tice, it assures a high quality weld 
from one side. Leaving an 0.050-in. 
land, at the inner surface a ;;-in. 
radius was machined and the joint 
face turned at 20 deg. for a little over 
one inch. The remainder of the sur- 
face was cut at a 10-deg. taper. 

After visual inspection of the joint, 
the type 347 stainless steel EB insert 
was tack welded around the 0.050-in. 
land of one of the hemispheres. This 
insert, manufactured by Arcos Corp., 
Philadelphia, is in the form of wire 
with a special cross sectional shape of 
amushroom. The head of this mush- 
room shaped wire was toward the in- 
side of the sphere. 

The two pressure vessel halves were 
placed in a holding fixture that care- 
fully located the hemispheres to per- 
mit the tack welding of the insert to 
the other side of the joint. The 
homogeneous reactor sphere was 
purged and filled with an inert gas 
to protect the inner surface against 
oxidation during the root pass and 
the laying of the next two beads. 

‘All welding was done from the top 
by putting the holding fixture in the 
vertical position. Fusing of the weld 
insert by inert gas shielded tungsten 
arc from one side only proceeded in 
12-in. segments, approximately one 
eighth of the circumference. One 
segment was fused and the second 
bead was overlayed before going to 
the next segment. The third bead 
was run over both sections before the 
fixture was rotated for welding of the 
two diametrically opposite segments. 

After completion of the root pass 
and the two subsequent passes, the 
closure was subjected to radiographic 
inspection with no defects of any 
sort permitted. After the inspection 
proved satisfactory, subsequent passes 
were made with type 347 electrodes to 
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join the stainless steel to within }-in. 
of the interface of the stainless steel 
cladding and the mild steel where type 
308L was used for five passes. Armco 
iron was used for 12 passes and the 
remainder of the joint was completed 
with E-7016 electrodes. 

Key to the entire weld, performed 
from one side, is the root pass made 
with the EB weld insert. This tech- 
nique produced a smooth, defect-free 
base obtainable by no other process. 


Aluminized Suit to Protect Fire- 
Fighters.—An expendable aluminized 
suit for the protection of fire fighters 
and a unique emergency water supply 
system were used successfully in large 
scale fire tests conducted recently in 
northern California forests by the 
U. S. Army Research and Develop- 
ment Laboratories, Fort Belvoir, Va. 

Preliminary reports on the tests, 
held to evaluate fire fighting equip- 
ment for troop use, indicate that men 
wearing the suit over duty uniforms 
were able to stand within two feet of 
a forest fire for two to three minutes 
without discomfort, although the heat 
was so intense that exposed portions 
of helmet liners were blistered. 

Designed by the Engineer Labora- 
tories and fabricated by the Quarter- 
master Corps, the expendable suit is 
made of flame-retardant treated alu- 
minized kraft paper. It consists of a 
parka with hood, a face mask, a pair 
of leg sleeves, and a pair of mittens. 

Other results of the test conducted 
in cooperation with troops from Fort 
Lewis, Washington, showed the feasi- 
bility of an emergency water supply 
system, also designed by the Labora- 
tories’ Fire Fighting Branch. This 
consists of lightweight piping and 
hose, and an electrically-driven sub- 
mersible pump. 

Soldiers inexperienced in fire-fight- 
ing were able to couple 1000 ft. of pipe 
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in three and a quarter minutes and 
the lightweight hose was laid success- 
fully by a helicopter flying at tree-top 
level. 


Helium Wind Tunnel.—Helium, 


the gas that fillsslow-moving blimps, is 
being used in studies of problems 
associated with the ultra-high-speed 
flight of satellites, missiles, and mete- 
orites at the General Electric Research 
Laboratory. Dr. Guy Suits, GE vice 
president and director of research, has 
disclosed that the Research Labora- 
tory is now operating what is believed 
to be the largest and fastest helium 
wind tunnel in the world. 

The new facility, in which helium 
serves as a “stand-in” for air, has 
produced the highest Mach numbers 
ever attained in a wind tunnel, reach- 
ing velocities in helium corresponding 
to Mach number 28 (28 times the 
speed of sound), Suits revealed. The 
attainment of these record-breaking 
Mach numbers, he added, represents 
“progress toward better understand- 
ing of the problems of high-speed 
flight.”’ 

Speeds corresponding to a Mach 
number of approximately 25 in air 
were needed to launch the first U. S. 
satellite, ‘Explorer A Mach 
number of 35, based on sea level sound 
velocity, would be necessary for an 
object to escape the earth’s gravita- 
tional field. Mach numbers in this 
region may be reached in future 
experiments. 

The studies at the GE Research 
Laboratory, partly sponsored by the 
U. S. Air Force, make use of a metal 
tube in which helium flows past a 
model representing, for example, a 
missile or its nose cone. In most con- 
ventional wind tunnels, huge com- 
pressors send a stream of air past the 
model being studied. In this helium 
tunnel, the compressed gas is stored 


in tanks and then allowed to escape 
through a nozzle into the test section, 
in which a near-vacuum is maintained. 
As the expanding gas rushes down the 
tube past the model, very high veloc- 
ities, together with very low tempera- 
tures and pressures, are produced. 
Static temperatures only two or three 
degrees above absolute zero have been 
obtained at the Research Laboratory 
during runs lasting from seconds to 
minutes. 

“The growing interest of the Gen- 
eral Electric Company in very high 
speed flight has resulted in a pressing 
need for a greater understanding of 
the flow of air above Mach number 
20,” according to Robert H. Johnson, 
mechanical engineer conducting the 
helium tunnel experiments in the lab- 
oratory’s recently completed $1.2 
million Gas Dynamics Facility. 

“This interest in the mechanics of 
fluids is not, however, a recent devel- 
opment for General Electric,’’ Johnson 
added. “It began with the cut-and- 
try methods of the early days of the 
steam turbine and the oscillating wall 
fan. Today this empirical approach 
has matured into sophisticated theo- 
ries that make possible continued 
improvements in the design of steam 
and gas turbines, as well as the acqui- 
sition of new knowledge in the field 
of very-high-speed flight.”’ 

Helium is used in the tunnel in 
order to avoid some of the complexi- 
ties that face the researcher studying 
high-speed flight in air. At extremely 
high speeds, air breaks down into an 
assortment of elements, compounds, 
and ionized particles that bears almost 
no resemblance to the air we breathe. 
Helium, however, the stablest sub- 
stance in the universe, is not subject 
to the great changes that occur in 
air’s properties and composition. In 
addition, helium has an extremely low 
liquefaction temperature, which en- 
sures that it will not liquefy (as air 
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would) as it passes through the tunnel 
passages and through the shock waves 
that occur in front of the model. 

As a result of helium’s more tracta- 
ble behavior, it is possible, by analogy, 
to obtain information on the dynamic 
and viscous behavior of air at very 
high Mach numbers. This informa- 
tion, together with data from comple- 
mentary research tools (such as shock 
tubes), contributes to a better under- 
standing of the complex actual be- 
havior of air under similar conditions. 

Test sections with diameters of 
either 3.6 or 6 in. can be used in the 
tunnel, depending upon the Mach 
number desired. Models up to one 
inch in diameter can be studied, with 
data taken by means of schlieren 
photographs and manometer board 
readings. 


Circulation of Atlantic Studied.— 
Scientists at the Woods Hole (Mass.) 
Oceanographic Institution are analyz- 
ing extensive data on deep currents of 
the Caribbean Sea and the Atlantic 
Ocean brought back from a two 
months’ cruise by the research vessel 
Crawford. 

The information collected on the 
Crawford, when considered with other 
Atlantic 1GY data, will shed impor- 
tant light on the circulation of the 
Atlantic Ocean and particularly on 
the “age’’ of deep Caribbean water, 
according to Dr. Columbus O’D 
Iselin, Woods Hole director. The 
Woods Hole effort is part of the Na- 
tional Academy of Sciences’ Interna- 
tional Geophysical Year Committee 
program to enlarge man’s knowledge 
of the earth. 

The six-man scientific party on 
board the ship, headed by oceanogra- 
pher William G. Metcalf, obtained 
more than 2000 deep temperature 
measurements and an equal number 
of water samples which were analyzed 
for salinity, oxygen content and other 
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chemical constituents. The samples 
were obtained from the surface down 
to depths as great as 23,000 feet while 
the ship was stopped. In addition, 
some 357 temperature measurements 
from the surface to a depth of 900 feet 
were obtained while the ship was 
under way. 

Detailed examination of the obser- 
vations, and comparison with similar 
data obtained 30 years ago by the 
Institution’s R. V. Adlantis, will help 
to reveal how old the water is in the 
enclosed Caribbean Basin. Impor- 
tant are the depths of the “‘sills’’ in the 
passages between Caribbean islands, 
through which deep Atlantic water 
cannot enter under present-day cli- 
matic conditions. Sills are submarine 
ridges connecting two land masses; 
sill depth is the distance from the 
surface of the water to the top of the 
ridge. 

Using the recentiy-developed Woods 
Hole precision echosounder recorder, 
the Crawford determined that the 
Windward Passage between Cuba and 
Haiti has a sill depth of 5,046 feet. 
This passage and the Anegada Pas- 
sage in the Virgin Islands are the 
deepest sills leading to the Caribbean 
and prevent exchange between the 
coldest Atlantic water and the deep 
Caribbean water. 

The echosounder was constantly in 
use and obtained accurate bottom 
profiles along the entire track of the 
Crawford, ranging from Bermuda to 
the Virgin Islands to Venezuela, from 
Colombia to Haiti, Jamaica to Cen- 
tral America, back to Cuba and thence 
to Panama. 

To date in the Woods Hole-IGY 
Deep Current Program, more than 
7000 deep temperature and water 
samples have been taken as well 
as some 3760 shallow temperatures 
(down to 900 feet) and bottom pro- 
files along all tracks of the ships. 
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Through the combined efforts of the 
(British) National Institute of Ocean- 
ography and the Woods Hole Ocean- 
ographic Institution, a countercurrent 
under the Gulf Stream was charted 
last year. Eight crossings of the 
North and Central Atlantic Ocean, 
with some 140 “deep stations’ com- 
pleted en route, also were made by 
ships of the two institutions. The 
data from these cruises will shortly be 
forwarded to the IGY Data Centers. 


Helium Separation and Purifica- 
tion.—A diffusion technique which 
promises to facilitate greatly the large- 
scale separation of helium from gas- 
eous mixtures such as natural gas has 
been developed by K. B. McAfee of 
Bell Telephone Laboratories. He- 


lium with no detectable impurities 
may also be obtained by using this 
technique as a purification process. 


The demand for helium in the U. S. 
has risen much faster than the ability 
to produce it. Our only sources of 
supply are a few natural gas wells, 
where helium exists at concentrations 
ranging from less than 1 to 6 or 7 per 
cent. Most of this helium at present 
is wasted when the natural gas is 
burned for fuel. Any process for 
readily recovering helium would be of 
tremendous commercial value. 

In the newly developed diffusion 
process, helium is separated from a 
gaseous mixture, such as natural gas, 
by passing the mixture over the sur- 
face of glass which has a high permea- 
bility to helium and low permeability 
to other gases. Silica glass has been 
found to be particularly advantageous 
in this respect, having a permeability 
to helium as much as one thousand 
times greater than for hydrogen, the 
next most diffusible material. The 
permeability is high enough to permit 
relatively large quantities of helium 
to diffuse through the glass under 
proper conditions. 


To obtain appreciable quantities of 
helium, a large surface of glass must 
be exposed to the mixture, the glass 
must be very thin, and a high pressure 
differential should be maintained be- 
tween the two sides of the glass. An 
excellent configuration for providing 
such conditions consists of a bundle 
of fine capillary glass tubing so ar- 
ranged that the gas mixture flows 
around the outside of the tubing under 
high pressure and the helium is re- 
covered from the inside of the cap- 
illaries. Appropriate seals have been 
developed to facilitate such an 
arrangement. 

Silica or pyrex tubing can be drawn 
having an external diameter of two 
mils and a wall thickness of two- 
tenths of a mil. Such tubing can 
withstand a compressive stress in 
excess of one thousand atmospheres, 
and is ideally suited for the separation 
or purification of helium. Diffusion 
through the glass increases rapidly 
with temperature, and the tubing will 
withstand temperatures of 400° C. 
and higher over long periods of time 
without deteriorating. 

Tests on an experimental thin wall 
capillary diffusion cell made up of a 
bundle of capillaries indicate that a 
similar cell containing enough capil- 
laries to occupy about two cubic yards 
would pass nearly 1000 cubic feet per 
day of helium at room temperature 
with a pressure differential of 1000 
atmospheres, assuming a concentra- 
tion of 1 per cent helium in the gas 
mixture can be maintained. Increas- 
ing the temperature to 400° C. would 
permit the recovery of 100,000 cu. ft. 
of helium per day. A cell of this type 
will not deteriorate with time, and 
might be placed directly in a gas pipe- 
line. Thus it appears that this proc- 
ess may have a truly large scale 
potential. 

Exceptionally pure helium results 
from a single diffusion step, even 
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though the initial mixture contains 
as much as 90 per cent hydrogen and 
only 10 per cent helium. A careful 
analysis indicates that there is less 
than 0.0009 per cent hydrogen in the 
helium purified by such a step. He- 
lium purified by diffusion through 
glass by this process is among the 
purest gases which it is possible to 
produce. 

Results already obtained indicate 
that this process can be adapted to 
large scale commercial purification of 
helium, and to certain high volume 
commercial separations involving this 
gas. 


International Lighting Vocabulary. 
—aA new international vocabulary of 
lighting terms, the culmination of 20 
years of study by a working party 
of the Commission Internationale 
de l’Eclairage (C.I.E.), is now ready 
for distribution through the organiza- 
tion’s United States National Com- 
mittee. Printed in three languages— 
French, English, and German—the 
vocabulary was created for the ex- 
press purpose of facilitating commu- 
nication among scientists the world 
over who are working in the field of 
illumination. 

The publication, entitled “‘Interna- 
tional Lighting Vocabulary of the 
International Commission on I]lumi- 
nation, Volume 1 (2nd Edition),” 
contains 530 terms, with definitions, 
as well as numerous symbols and for- 
mulae. Its principal sections cover 
Radiation, Photometry, Colorimetry, 
Eye and Vision, Production of Light, 
Lamps, Auxiliary Apparatus, Lighting 
Techniques, Lighting Fittings, and 
Light Signals. 

Volume 2 of the publication will be 
published in 1959. It will contain 
the same terms, without definitions, 
and will be printed in ten languages. 
These are: French, English, German, 
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Danish, Dutch, Italian, Polish, Rus- 
sian, Spanish, and Swedish. The ma- 
terial will be arranged systematically, 
with ten alphabetical indexes. 

The United States National Com- 
mittee has placed a prepublication 
order for the first volume and is offer- 
ing it at a reduced price. A.limited 
number of copies are now available 
at $2.50. It is expected that the 
price will be raised when the initial 
supply is exhausted. 

Orders for the publication with re- 
mittances payable to “‘U. S. National 
Committee, C.I.E.”’ should be sent to 
Mr. T. D. Wakefield, Treasurer, 
U.S.N.C., Vermilion, Ohio. 


Laboratory Automation System.— 
Complete automation in the chemical 
laboratory will be introduced in the 
United States this year. It will be 
made possible by a system of instru- 
ments that will perform the most 
exacting scientific tests, record their 
findings, and take proper corrective 
action when necessary, all without 
human help. The system will be 
offered this fall by Chicago Apparatus 
Company. 

Called the Analmatic, the instru- 
ments are made by Baird & Tatlock 
Ltd., Great Britain, and are exten- 
sions of a basic system the company 
designed and built for the British 
Atomic Energy Authority. The sys- 
tem completely frees the scientist from 
routine analytical work. 

The Analmatic will be displayed for 
the first time in this country at the 
National Chemical Exposition by the 
Chicago section of ACS Sept. 9-12 in 
Chicago’s International Amphitheatre. 
Chicago Apparatus has been named 
exclusive United States and Canadian 
distributor for the Analmatic. 

The system is far more comprehen- 
sive automation than anything that 
has been tried in American labora- 
tories. The term “automation” very 
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often is used too loosely to describe 
single instruments that give a quick 
reading or short cut the analysis pro- 
cedure, once the scientist has prepared 
the testing sample. 

Here are the operations performed 
by the Analmatic instrument system 
without human help: 


1. Collects a succession of exactly 
measured samples of the material to 
be analyzed. 

2. Mixes the successive samples 
with metered amounts of reagents, as 
necessary. For example, analyzing 
uranium concentrates involves the ad- 
dition of ammonium thiocyanate. 

3. Heats or agitates the solution as 
required. 

4. Performs the necessary test by 
color measurement, by titration to 
pH, by potentiometric or ‘dead stop” 
end points, by photometric measure- 
ment, etc. 

5. Records the test results as graph 
or dot diagrams, or in digital form. 

6. Thoroughly washes all the pi- 
pettes and other vessels used to hold 
liquids, insuring that there is no con- 
tamination that will affect subsequent 
tests. 

7. If the analysis shows that the 
test material deviates from a given 
critical value or range, a visual or 
aural alarm is given. If desired, the 
alarm can be replaced by a control 
system that will adjust the production 
process automatically, compensating 
for the previous error. 


The instrument system, not an as- 
sembly line item, is built with typical 
British craftsmanship. Price of the 
instruments will depend on individual 
requirements, since each system is en- 
gineered for the user. Probable price 
range will be from $15,000 up to 
$50,000. The reason the range may 
be so wide is that the system will vary 
from the relatively simple to the very 
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complicated. In a given case, Baird 
& Tatlock can engineer a system cap- 
able of handling several hundred daily 
determinations of a complicated chem- 
ical assay. In other instances, a rela- 
tively simple and less expensive Anal- 
matic is now used by several English 
firms for the rapid and accurate de- 
termination of moisture by the Karl 
Fischer reagent test. 

The first analyzer was put in opera- 
tion by the British Atomic Energy 
Authority in South Africa to check the 
waste flow in a uranium purifying 
process, guarding against too much 
uranium being “washed down the 
drain.” In this application, the Anal- 
matic does the work of two skilled 
analysts who were taking samples at 
15-min. intervals. For 24-hr. opera- 
tion, the manpower saving is tripled. 
Other current users in England include 
a manufacturer of penicillin and other 
pharmaceuticals, the largest chemical 
firm in the country, and a petroleum 
products plant. 

According to Chicago Apparatus, 
the following industries in this country 
are the most likely to have automatic 
laboratory applications: pharmaceu- 
ticals, petroleum processing, many 
different kinds of chemical production, 
paper manufacturing, food processing, 
plastics production. 

Analmatic equipment is at present 
available only for clear liquids. Baird 
& Tatlock has designed apparatus for 
external filtration of samples and 
reagents entering the analyzer, but 
has not worked on automation tech- 
niques for testing slurries, suspensions, 
or solids. 

As American agent, Chicago Ap- 
paratus will provide technical assis- 
tance for all installations. Actual 
engineering systems for the first few 
months will be performed by Baird & 
Tatlock personnel who will come to 
this country in August. 
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The Franklin 
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Development 


One of the nation’s largest 
industrial research laboratories 


Where a fully qualified staff of engineers 
and scientists works as a competent and 
versatile research team—bringing a fresh 
scientific approach to problems of materials, 
methods and machines of modern industry. 


in the fields of 


Chemistry and Physics * Electrical Engineering 
Mechanical Engineering . Solid State Physics 
Nuclear Engineering 
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